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Abstract 

We provide an extension of a recent approach to study non-equilibrium thermodynamics [Phys. 
Rev. E 81, 051130 (2010), to be denoted by I in this work] to inhomogeneous systems by con- 
sidering the latter to be composed of quasi-independent subsystems. The system S along with 
the (macroscopically extremely large) medium S forms an isolated system Eq. Starting from the 
Gibbsian formulation of the entropy for Sq, which is valid even when Sq is out of equilibrium, 
we derive the Gibbsian formulation of the entropy of E, which need not be in equilibrium. We 
show that the additivity of entropy requires quasi-independence of the subsystems, which requires 
that the interaction energies between different subsystems must be negligible so that the energy 
also becomes additive. The thermodynamic potentials of subsystems such as the Gibbs free energy 
that continuously decrease during approach to equilibrium are determined by the field parameters 
(temperature, pressure, etc.) of the medium and exist no matter how far the subsystems are out 
of equilibrium so that their field variables may not even exist. This and the requirement of quasi- 
independence make our approach different from the conventional approach due to de Groot and 
others as discussed in the text. As the energy depends on the frame of reference, the thermody- 
namic potentials and Gibbs fundamental relation, but not the entropy, depend on the frame of 
reference. The possibility of relative motion between subsystems described by their net linear and 
angular momenta gives rise to viscous dissipation. The concept of internal equilibrium introduced 
in I is developed further here and its important consequences are discussed for inhomogeneous 
systems. The concept of internal variables (various examples are given in the text) as variables 
that cannot be controlled by the observer for non-equilibrium evolution is also discussed. They are 
important because the internal equilibrium in the presence of internal variables is lost if internal 
variables are not used in thermodynamics. It is argued that their affinity vanishes only in equilib- 
rium. Gibbs fundamental relation, thermodynamic potentials and irreversible entropy generation 
are expressed in terms of observables and internal variables. We use these relations to eventually 
formulate the non-equilibrium thermodynamics of inhomogeneous systems. We also briefly discuss 
the case when bodies form an isolated system without any medium to obtain their irreversible 
contributions and show that this case is no different than when bodies are in an extremely large 
medium. 
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I. INTRODUCTION 



A. Nature of the Problem 

In an earlier paper which we will refer to as I in this work, we have considered some 
of the consequences of applying the second law of thermodynamics to an isolated system Eq, 
which consists of a macroscopic system of interest S containing a fixed number of particles 
(atoms or molecules) surrounded by an extremely large medium S; see Fig. [TJ From now 
on, it will be implicitly assumed that the medium is extremely large to be unaffected by 
the system. (Later in this work, we will also consider S to have a fixed volume V instead 
of N or fixed V and A^. Furthermore, we will also consider briefly the case of many similar 
size systems forming the isolated system Eq without the extremely large medium S.) The 
motivation has been to develop this approach to obtain a non-equilibrium thermodynamic 
description of the open system E under various conditions. Throughout this work, we will 
use body to refer to any of the three kinds systems: the isolated system, the medium or 
the system. All quantities related to Sq, S and S will be denoted by a suffix 0,~ over the 
top, and without any suffix, respectively. All quantities related to a body will be denoted 
without any suffix, as if we are dealing with an open system. Similarly, in this work, we will 
say that the system is open when it is in a medium. Even though it is not the common 
usage, this should not cause any confusion as the context will be clear. 

To avoid complications due to external shear, we had only considered S under no external 
shear in I. This restriction is easily removed as we will do here. The isolated system will 
still have no external force acting on it to ensure its isolation; see below also. We will now 
allow forces acting at the surface dV of the system or any of its subsystems; see Fig. [2](a). 
These forces must balance the internal stress tensor for mechanical equilibrium. Thus, the 
force tidf acting on a surface element df must equal the stress force TijUjdf (summation over 
repeated indices implied), and we have for the surface traction force 

ti = TijUj. (1) 

Here, n is the outward unit normal at the surface element df surrounding a point on the 
surface. This condition must be satisfied at every point on the surface dV. The net force 
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FIG. 1: Schematic representation of a macroscopically large system S and the medium S sur- 
rounding it to form an isolated system Sq. The system is a very small part of Sq. The medium 
is described by its fixed fields To,Po, etc. while the system, if in internal equilibrium (see text) is 
characterized by T{t),P{t), etc. 



and torque acting on the system are given by 

F = ^tdf, K = y"r(") X tdf, 

dV dV 

respectively; here, for convenience, r*^^^ is taken to be the radius vector of the surface element 
with respect to the center of mass of the system. The external forces are responsible for the 
deformation of the system, and also result in the translation and rotation of the system. Let 
us consider an infinitesimal volume element dV of mass dm, which is moving with a velocity 
v(t) and rotating with an angular velocity ijj{t) and has an intrinsic angular momentum 
m{t)dV . The linear and angular momenta of the system in some fixed frame are given by 

MY{t) =Jy^^^v{t)dm, M(t)=MR(t) x V{t)+Jy^^^r{t) x v{t)dm + Jy^^^m{t)dV, 

respectively; here M = rriQN is the mass of the system ( mo being the mass of a particle), 
which we consider fixed for fixed A^, and R(t) and V(t) are the location and the translational 
velocity of the center of mass in this frame. If the frame is taken to be the center of mass 
frame, then R(t) and V(t) are zero. If the body as a whole is stationary, then the energy 
of the body is known as the internal energy. Such a stationary situation was considered in 
our previous work where no motions were considered. This limitation will be removed 
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FIG. 2: We show schemtically the two subsystems o"i and 02 (72 (t) > Ti[t)) forming the system 
S in an extensively large medium in (a) and by themselves forming an isolated system without 
an extensively large medium in (b). The heat output dQ{t) in (a) by 02 is the sum of dQ'{t) 
and dQ2{t), while the heat intake by cJi is the sum of dQ'{t) and dQi{t). As we are dealing with 
an isolated system, the heat input and output must be equal. Therefore, dQi{t) = dQ2{t). The 
equality of the heat input and output is also true in (b). As the heat transfer occurs between objects 
does not occur isothermally, there is irreversible entropy generation due to each heat transfer. We 
will study this issue later in Sect. IXI B[ 



here so as to allow for relative motions (translation and rotation) between the system and 
the medium. These relative motions are the additional sources of viscous dissipation and 
give rise to additional irreversibility. The irreversibility due to temperature difference (such 
as between the system and the medium) has already been considered in I. 

We should mention here the recent somewhat comprehensive investigation carried out by 
Bouchbinder and Langer J] who also consider a system under external shear; however, our 
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approaches and emphases are very different. We should also mention earlier very different 
equilibrium-like attempts by Lubchenko and Wolynes jsl and by Ottinger ^\. Mention 
should also be made of a very interesting phenomenological approach by Oono and Paniconi 
0] on steady state thermodynamics, which was later advanced by Sasa and Tasaki [§1. The 



classical local non-equilibrium thermodynamics due to de Bonder |9Nl2| is close in spirit to 



our approach and that of Bouchbinder and Langer W. It will be the standard formulation 
with which we will compare and contrast our approach initiated in Therefore, for 
the sake of continuity and clarity, we briefly discuss the classical formulation involving 
local equilibrium in Sect. HIl We note that there are other versions of non-equilibrium 
thermodynamics usually known as the extended, rational and GENERIC non-equilibrium 



thermodynamics 
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14| : however, we do not discuss these formulations in this work. 
It is well known that the second law of thermodynamics for the isolated system Eq states 
that 

dSr)(t) , . 

^ > 0, (2) 

where 5*0 (t) denotes the entropy of Sq at some instant t. It should be stressed that the 
isolation of Sq requires that there be something outside of Sq from which it is isolated. 
Therefore, we will assumes that Sq is confined to a finite though extremely large volume Vq 
[isl l . The isolation requires that we neglect all interactions, such as gravitational interactions, 
of Sq with the outside that cannot be shielded. All interactions with the outside should 
be relatively very weak to be negligible. All relevant interactions must be confined inside 
the volume Vq. This to not to be taken as a weakness of our approach as including the 
interactions with outside will only make Sq an open system, so that our investigation of an 
open system, which is our primary concern, can then be applied to it. 

We should emphasize the following important point. The law of increase of entropy in 
Eq. ([2]) should give a pause to those readers who believe that the concept of entropy is 
meaningful only for an equilibrium state and that the entropy cannot be defined for a non- 
equilibrium state. This belief is unfounded. The mere fact that the second law dictates the 
approach to equilibrium clearly implies that the entropy exists even when the system is out 
of equilibrium. This issue and the related history, in particular the Gibbs formulation of the 
non-equilibrium entropy, has been reviewed recently [isj, and references cited therein], and 
we refer the reader to this for further details. We should stress that the Gibbs formulation 
of the entropy requires that the dynamics of the system and that of the isolated system are 



not deterministic, but rather stochastic; see 



15| for further details. 



B. Important Restrictions in I and Their Removal in the Current Work 

An assumption that was imphcit, but not stated, in I was that there was no relative 
motion between the system and the medium and that the isolated system was stationary. 
This is normally the case in practice in which the various bodies (Sq, S and S) are stationary 
in the laboratory frame of reference; the latter we will denote by L in this work. In general, 
one can allow for a translation and rotation of a body, as done here because we now wish to 
study their effects on the deformation of the system. Because of the stationary assumption, 
there is no difference between the energy and the internal energy for Sq, S and S. This 
is a very common but useful assumption as the entropies depend on the internal energies 



and not on the energies [16|. The latter energies may contain the contribution from the 
translation and rotation of the system as a whole. In contrast, the internal energies are the 
energies systems have in a frame in which they are stationary. Whenever we discuss both 
energies together in the following, the internal energy will be denoted by a superscript "i" 
to distinguish it from the energy, which is denoted without the superscript; otherwise, it 
will be clear which energy we are considering. For a stationary system in the lab frame £, 
the internal energy is the same as the energy. The translation of a body as a whole merely 
affects the energy, but not its thermodynamic properties. However, the rotation of a body 
as a whole gives rise to centrifugal potential energy that modifies the energies of microstates 



and has to be carefully incorporated in any thermodynamic investigation 
as done in Sect. IVIIIj see also Appendices lA] and [Bl 
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see Sect. 34] 



C. Present Goal and the Layout 

We had focused exclusively on the system as a whole without worrying about its internal 
parts in I. As entropy palys a central role in our development, we discuss the formulation 
of entropy and thermodynamic averages in Sect. IIIIl where we show the entropy itself as a 
thermodynamic average. Under the assumption of internal equilibrium in the system, there 
were no relative motions between its various parts that could add additional irreversible 
processes inside the system. In this work, we will remove this limitation and treat the system 
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E as inhomogeneous as glasses normally exhibit both spatial and temporal inhomogeneity; 



see 



17H19| for some recent investigations. We do this by considering S to be composed of a 



collection of a large number A^^g of subsystems ak, k = 1, 2, ■ ■ ■ ,Ns, which may be different 
from each other to allow for inhomogeneity and for relative motions and shear forces between 
different subsystems in terms of internal variables. Each subsystem is still macro scopically 
large so that we can introduce a legitimate entropy function s^; see Sect. IX Bl for further 
elaboration. 

Apart from the observables that can be manipulated by the observer, there also appear 
internal variables often used in describing glasses, as is well known from the early works 
of Davies and Jones 
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2l| . The latter variables cannot be manipulated by the observer 



and were briefly introduced in I, but not explored. We will also remedy this situation here 
and consider glass as an inhomogeneous system with internal variables. Specifically, we 
treat translations and rotations of various parts of a system as internal variables that are 
generated by surface traction forces. The alternative approach is to use the traction forces 
and the strains instead; see for example ^. The phenomenological ideology introduced by 



Davies and Jones 



201], which has been recently reviewed by Ottinger |6i], is by now standard 



and has been discussed in several textbooks; see for example [22|, |23j. The observables and 
internal variables will be collectively called state variables; see Sect. Illll for proper definitions 
of these terms. 

In Sect. ??, we discuss the consequences of internal equilibrium and its similarity with and 
differences from the concept of local equilibrium js jj] discussed in Sect. [Tll In particular, 
we argue in the form of Theorem [T] that the system can only sustain a uniform translation and 
rotation in internal equilibrium. It is assumed here that there are no additional conditions 
(such as the potential flow in a superfiuid) on the velocity. The proof is trivial but the 
theorem has far reaching consequences. For example, it implies that the uniform rotation 
must be about a principle axis of inertia. A simple way to understand internal equilibrium 
is to think of the system as follows. We first disconnect it form the medium with which 
it is not in equilibrium, and connect it to another medium whose conjugate field variables 
(T, P, etc.) are exactly the same as that of the system. The system remains in equilibrium 
with this medium so that there will be no irreversible process in the system. We discuss 
the generalization of equilibrium Maxwell relation to systems in internal equilibrium. The 
condition for additivity and quasi- independence is considered in Sect. El where we prove 
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that the Gibbs entropy formulation in Eq. (12T]) is also applicable to an open system, which 
is assumed to be quasi-independent of the medium. Various thermodynamic potentials 
are identified in Sect. |Vl] that are in accordance with the second law. It is here that we 
see a clear distinction between our approach and the conventional non-equilibrium theory 
exploiting the local equilibrium concept js iSl. Internal variables are discussed in Sect. IVIII 



We prove that the chemical potential or the affinity associated with an internal variable must 
be zero when the system is in equilibrium. We also prove that the entropy expressed solely in 
terms of observables when there are independent internal variables must explicitly depend on 
time so that while the system is in internal equilibrium with respect to all state variables, 
it is not in internal equilibrium with respect to only observables. A system undergoing 
uniform translation and rotation is studied in Sect. IVIIII where we also develop the Gibbs 
fundamental relation for such a system. We then apply the results from this study to an 
isolated system in which the system and the medium undergo relative translational motion in 
Sect. IIXI but the system is homogeneous. An inhomogeneous system with relative motions 
between its subsystems is studied in the next section. We also discuss in this section the 
case of several bodies, each in internal equilibrium but different from others, that form an 
isolated body together by themselves without a medium; see Fig. [2](b). These bodies are 
macroscopic in size, but are not extensively as large as a medium. We discuss a direct method 
of calculating the irreversible entropy generation in each body in terms of the equilibrium 
state of all the bodies. We find that the same results are also obtained by bringing all bodies 
in contact with a medium as shown in Fig. EJ^a). This equivalence is used to prove Theorem 
[21 All these investigations are extended to include extra observables and internal variables. 
The final section contains concluding discussion and a brief summary of our results. 

II. LOCAL NON-EQUILIBRIUM THERMODYNAMICS: A BRIEF REVIEW 

As entropy is the central quantity appearing in the second law in Eq. ([2]), we will pay 
close attention to its determination, although this is usually not done in classical local 



non-equilibrium thermodynamics [9l-ll2l|. where ^ts existence is taken as a postulate with an 
implicit assumption that it is always additive [9]. The entropy 5 of a system is defined in 
terms of the local entropy density s(r) per unit volume: 

S ^ jys{r)dV. (3) 
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The local temperature T and pressure P are assumed continuous functions of the location 
r and time t, and are postulated to always exist. The additivity of the energy E results in 

1 



e(r) + -mop(r)v (r) + m(r)-a;(r) + ■(/'(r) 



dV. (4) 



Here, e{r)dV is the internal energy, |mop(r)v^(r)(iy and m{r)-uj{r)dV the translational and 
rotational kinetic energy density, respectively, and ip{r)dV the additional energy contribution 
due to interactions not included in the internal energy density in a volume dV of the system; 
the local mass and angular momentum densities are given by mop and m(r), respectively; 
compare with Eq. flll2p derived later after limiting it to the volume element dV. 

The functional form of the entropy density depends on the nature of the system. For 
example, for a simple system containing a fixed number of structureless particles, it is 
assumed to be a function only of the internal energy density e and the local number density 
p; see for example, [9|, see Eq. (III. 14)] 

s(r) = .(e(r),p(r)). (5) 

The local Gibbs free energy density ^ is given by 

g = e-Ts + P. (6) 

whether local equilibrium exists or not. (The unconventional use of the symbol ^ instead of g 
will become clear later when we discuss the Gibbs free energy, which follows from the second 
law and which continuously decrease as the system approaches equilibrium.) However, no 
direct method of calculating the entropy is given in this approach except by assuming as 
another postulate the validity of Gibbs fundamental relation, which for a simple system with 
no internal variables reads 

Td{s/p) = d{e/p) + Pd{l/py, (7) 

this postulate is a consequence of assuming local equilibrium j^, see Eq. (HI. 15)]. The 
presence of the local temperature T and pressure P in the fundamental relation imposes 
strong conditions on the nature of the entropy in that its partial derivatives are related to 
the given T and P under local equilibrium, which follow from Eq. ([71). 



We avoid the above issues in the conventional non-equilibrium thermodynamics 9l-ll2| 
by first identifying the entropy of the system in terms of microstate probabilities as de- 
scribed below, see Eq. (|2T1) . and then use the concept of internal equilibrium to introduce 
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the temperature and pressure in terms of this entropy; the latter are defined only when 
there is internal equilibrium This thus avoids the issues inherent in the conventional 
approach. Our approach is not local in that we always deal with quantities S, E, etc. related 
to macroscopically large systems or subsystems as opposed to the conventional thermody- 
namics which deals with local quantities s, e, etc. As a consequence, not only S, E, etc. but 
also the temperatures, pressures, etc. associated with these systems or subsystems will not 
always be continuous functions of space at the interfaces. Thus, we will not impose continu- 
ity in space on these quantities, which makes our approach distinct from the traditional local 



non-equilibrium approach of de Bonder |9l-|l2| where these quantities are always treated as 



continuous. In the latter approach, a system can be broken into subsystems, each sufficiently 
small, yet large enough to be in internal equilibrium to satisfy Gibbs fundamental relation. 
This hypothesis is known as the local equilibrium hypothesis. We add another requirement, 
that of quasi-independence of the subsystems in our approach, which we believe to be ex- 
tremely important. Only this requirement ensures that the entropy retains the additivity 
property and also remains a state variable, as we will discuss later in Sect. |Vl 

A. Helmholtz Theorem 

To accommodate relative motion, we need to allow for surface traction forces ti, which 
then give rise to internal forces causing the deformation of the system. These forces can 
be related to the induced stress tensor Tij as shown in Eq. ([T]), and will result in a motion 
of the system due to non-zero net force and torque acting on the system. It is well known 
that the local motions for a deformable system can be described as a combination of three 



25|: 



distinct types of motions 2J, 

(a) a pure translation 

(b) a pure strain, and 

(c) a pure rotation 

by expressing the instantaneous difference in the velocity 6v at two nearby points sepa- 
rated by a displacement vector 6r as 

Svi = ipijSxj + Qij6xj, (8) 
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where the symmetric tensor 



_ 1 , dvi dvj , 



represents the rate of strain tensor and the antisymmetric tensor 

^ _ 1 dvj dvi _ ^ 

is the vorticity tensor, and represents the axial vector cD = |9 x v associated with Qij. The 
second term QijSxj in Eq. represents the components of the vector u) x 6r. One should 
think of 5v as the relative velocity between two neighboring points separated by 6r. 

The first contribution represents a pure straining motion while the second contribution 
represents a rigid-body rotation. For example, a simple shearing motion in which plane 
layers of the system slide over each other can be treated as a combination of a pure strain 
(with no rate of volume change) and a rotation 25|. 



B. Stress Tensor 

The motion at the local level can also be studied directly by considering the stress tensor. 
The stress tensor is normally expressed as a sum of the non-dissipative and dissipative or 
viscous contributions 

Tij = 0-ij + cr[j, (9) 

in which the viscous contribution cr^'j is some function that depends on the velocity gradients, 
i.e., on 

dvi/dxj , d'^Vi/dxjdxk, etc. 

Thus, we can express it as 

a'ij = AijMfkl{dVm/dXn,d^Vm/dXndXp, ...), 

where Aijki does not depend on the velocity distribution and fki is a function of the velocity 
gradients. For example, in a linear approximation using only dvi/dxj, fki is taken to be 

fki = fodvk/dxi 

with a constant /o, which could be conveniently absorbed in Aijki- In this approximation, 
we see from Eq. ([8]) that cr-^ depends on both the rate of strain tensor and the vorticity u). 
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In general, we can partition cij' • into a symmetric and an antisymmetric part, the latter due 



to the presence of intrinsic rotation of the system 



9j and describing the role of the rotational 



viscosity. Thus, we can also partition Tjj into a symmetric and an antisymmetric part. We 
refer the reader to Chapter 12 in the monograph of de Groot and Mazur P]. Of course, 
there may be symmetry reasons such as the isotropy of the system that would forbid the 
dependence on vorticity, in which case there would be no antisymmetric part. It is possible 
to show 9| that the rate of change of the intrinsic angular momentum is determined solely 
by the antisymmetric part of r/^. Thus, the absence of this part will imply the conservation 
of the orbital and intrinsic angular momentum separately. 

It can now be shown Eqs. (XII. 18) and (XII. 18)] that the antisymmetric part 

contributes a term proportional to 

r ■ (uj—uj) 

to the rate of change of the internal energy e and to the entropy production. Here, r is the 



vector associated with rf- 



and uj represents the angular velocity of rotation of the system; cf. Eq. (|4]). Thus, this 
contribution vanishes for uniform rotation Cj=(jj as expected. For Cj^uj, there would be 
precession of the local volume element j26| about the direction of cj, so that the rotational 
viscosity would play an important role until uj becomes equal to uj. 



C. Energy Balance and the First Law 

The decomposition in Eq. ([9]) allows us to break the surface traction also in two terms 
related to the individual contribution: 

t = cr + cr', 

where the two new terms are defined similar to that in Eq. ([1]). The non- viscous contribution 
is sometimes expressed by taking out the isotropic pressure term as follows 

cr = — Pn + £, 
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assuming that the thermodynamic pressure can be defined. The rate at which the work is 
done on the system is given by the average (over all microstates) of the stress tensor at the 
surface dV of the system: 

^= ft,vfdf= fa,vfdf+ fa'^vfdf. (10) 

dV{t) dV{t) dV{t) 

where dfi and v^^'^ are the components of di{t) and the velocity v*^*^) of the surface element, 
respectively. When the pressure can be defined, this rate can be expressed as 



dV(t) dV(t) 

In this work, we will not consider latent heats and chemical or nuclear reactions within 

the body. In this case, the rate at which the heat is added to the system is obtained by 

considering the heat fiux through the surface dV and is given by 

dQ{t) 
dt ' 

where dQ{t) is the heat added to the system in time dt. We can then write down for the 
rate of change of the energy due to the dynamics in the system as 

dm ^ dm . dw'jt) 

dt dt dt ' ^ ' 

which is a restatement of the first law of thermodynamics; see also ^. 

It should be remarked that the balance equation f|T2|) is in an integral form for the entire 
system and should be contrasted with the differential form commonly stated in textbooks; 
see for example ol. The latter formulation is valid for infinitesimal volumes. Here, we are not 
interested in such a local description. Our main focus is to consider regions of the system 
that are macroscopically large enough so that they can be treated as quasi- independent. 
Under this condition, the entropy of the system can be approximated to a high accuracy by 
adding the entropies of the subsystems. We discuss this point further in Sect. IIX A[ 



D. Need for Internal Equilibrium 

If the system as a whole is stationary, then the average velocity V()f;) = 0. In this case, 
the energy of the system is the internal energy. We will consider this case below in this 
section for simplicity. 
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The first law statement in Eq. (fT2|) does not by itself allow us to determine the way 
the entropy of the system changes. This law is valid even if the system is not in internal 
equilibrium. In the absence of internal equilibrium, there is no way to determine the change 
in the entropy from this law. We need to relate dE{t) with the change dS{t) in the entropy 
to determine the latter. However, if we now assume that the system is in internal equilibrium 
(see I and Sect. ?? below), then the entropy no longer has an explicit t-dependence. In this 
case, we can write down the differential of the entropy S{E, V, N) of a monatomic system 
of neutral particles with fixed iV as [ij 

dS = JrrdE + ^dV; (13) 
T{t) T{t) ' ^ ' 

compare this equation with Eq. ([7]). We can now use Eq. (fT2|) in this equation to determine 
the rate of change of the entropy. However, as we discuss in Sect. ??, in this case there is 
no additional irreversible entropy production arising from viscous stress tensor. Thus, the 
last contribution in Eqs. ( ITOl) or ( ITTj) vanishes |27!, Sect. 7.4.2]. If there is no shearing at the 
surface, then the second term in Eq. (ITT]) also vanishes, and we obtain the standard relation 



this identification was used in I 



oLS(t) _ dQify 
dt ~ dt ' ^ > 



E. Reversible Entropy Change and Irreversible Entropy Production 

We wish to emphasize, as was done in the previous work see Eqs. (16) and (18) in 
particular] that this assumption of internal equilibrium does not mean that the irreversible 
entropy production diS{t) in the system is absent unless the system happens to be in equi- 
librium with the medium. To see this most easily, we recount from [l| that 

,(E)^(,) = ^ (15) 

is the change in the entropy due to heat transfer dQ(t) to the system that is in internal equi- 
librium and depends on its instantaneous temperature T(t). We have used the superscript 
E to indicate that the entropy change we are discussing is due to energy (heat) transfer for 
which the associated conjugate variable, see Eq. (1271) . is i/e = {OS / dE)v^N = ^/T{t) and 
takes the value ?/oe = l/^o in equilibrium. It follows from this that the reversible entropy 
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change dl S{t) in the system depends on its equilibrium temperature Tq, which is also the 
constant temperature of the medium, and is given by 

d^J^'Sit) = (16) 

regardless of the instantaneous temperature of the system. This results in [l| 

df^S{t) = dQ{t) - ^) = F^it)dQ{t) > (17) 

in all cases. Here, 

F^{t)^F[yit)]^yit)-yo, (18) 

with y = yo representing the equilibrium value of the conjugate field y associated with the 
observable X, represents the thermodynamic force for the flow of X. Since T{t) in Eq. (ITB]) 
undergoes a change dT oc dQ{t) due to the heat transfer, the heat transfer is not isothermal. 
Thus, it should not be surprising that there is an irreversible entropy generation as part of it. 
On the other hand, the determination of deS{t) in Eq. (fT6|) requires the heat transfer to be 
isothermal for the process to be reversible. What is surprising is that df.S{t) is determined 
not by the current state of the system at t, but its equilibrium state in the future so that 
the concept of causality is inapplicable [15|]- It is this particular aspect of deS(t) that is the 
cornerstone of the second law of thermodynamics: diS{t) > also depends on the future. 
We will make use of this observations later. 

The situation with other extensive variables like volume V is no different. As shown in 
|l|, the entropy change and the reversible entropy change due to a change dV are given by 

d'-''^S{t) = P^^dVit); rfW5(t) = p^dVit); (19) 

It follows from this that the irreversible entropy generation df^^ S{t) is given by an identical 
formulation as above for d^^^ S{t) 

d^Sit) = dV{t) - ^) = dV{t)Fy{t) > 0, (20) 

with Fy{t) given by Eq. f|T8l) with yy = P{t)/T{t), see also Eq. fl28l) . represents the 
thermodynamic force for the "flow" of volume V . 

It should be stressed that the validity of Eq. fll3p follows from the continuity of the 
entropy and the existence of the internal equilibrium. This is a general relation and is valid 
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for all systems whose macrostate can be described by the three variables E, V, and A^. On 
the other hand, the form of Eq. f lT^ depends on particular form of the processes that go on 
in the system. Thus, it is process-specific; recall that we have eliminated various processes 
in deriving this equation. We can incorporate the missing additional contributions in Eq. 
f ll4p by properly introducing internal variables to describe the microstates of the system and 
generalizing Eq. ( !T3|) to include other extensive observables. The issue of internal variables 
is taken up in Sect. IVIIt and the generalization of Eq. f lT3|) is taken up in Sect JVIIIl 



III. ENTROPY AND AVERAGES 



A. Isolated System 

The entropy of an isolated system such as Sq, whether in equilibrium or not, is given by 
the Gibbs formulation 

So{t) = -j:p^{t)\np^{t), (21) 

a 

where Pa(t) is the time- dependent probability of the a-th microstate of the isolated system 



So; see a recent review 15| of this formulation for an isolated system, regardless of how 
far it is from equilibrium. It is assumed that the dynamics of the system is stochastic and 
not deterministic, as the latter dynamics makes the above entropy a constant of motion in 
direct contradiction with the second law. In a deterministic, i.e. Hamiltonian dynamics, a 
microstate uniquely evolves into a microstate, while in a stochastic dynamics, a microstates 
evolves into several microstates with certain probabilities. The unique evolution is time- 
reversible, which results in the entropy being a constant of motion. In a stochastic dynamic. 



the entropy can only increase [15|. Since the system is isolated, we do not allow external 
forces acting on it; we of course neglect weak stochastic forces acting on it. Thus, any defor- 
mation if it occurs is due to internal forces. We have defined the entropy as a dimensionless 
quantity, which is equivalent to setting the Boltzmann constant k-Q equal to 1. The collection 
a. = {a} of these microstates along with their non-zero probabilities represents a macrostate 
A4q or simply of Sq. As equilibrium and non-equilibrium thermodynamics is an exper- 
imental science, a macrostate of any body is specified in terms of a set of some extensive 
observables X that can be manipulated by the observer. The same set of observables are also 
used to identify the microstates. Most often, the microstates are identified by the energy, 
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volume, and the number of particles that form the elements of X, because of their special 
role in thermodynamics. If there are other extensive mechanical quantities (quantities that 
are not purely thermodynamic in nature) such as the total linear and angular momenta, 
then the microstates are characterized by all these extensive quantities, collectively denote 
by X; see the discussion in Sect. IVIII for more details. Apart from these observables, a body 
can also be specified by a set I of internal variables [9-12| that have been found very useful 
in describing glassy dynamics and will be discussed later in Sect. IVIII Indeed, one needs 
internal variables to also explain the time-evolution of an isolated system towards equilib- 
rium since all its observables remain constant. Thus, the macrostate of the isolated system 
can use the internal variables for its specification. We will take these internal variables to 
be also extensive and call both of them as state variables and collectively denote them by Z. 
Taking internal variables as extensive allows us to deal all state variables on equal footing, 
so that generalization from observables to internal variables becomes almost trivial. 

Let us continue with the discussion of the isolated system. In general, microstate prob- 
abilities Pa{t) are functions of the state variables Zq along with t. As a consequence, the 
entropy S'o(Zo(t), t) is also going to be a function of Zo(t) and t. There are situations, when 
the entropy can also depend on some external parameters such as the angular velocity of 
the rotation of the frame of reference. These parameters need not necessarily be extensive 
just as t is not. The number of state variables are too limited for a complete microscopic 
description of the system, but is sufficient to describe the macroscopic conditions of a system. 

The observables remain fixed for the isolated system Sq. The internal variable Io(t), if 
present in the isolated system, is normally a function of time; its time-dependence describes 
the temporal evolution of Sq. While the microstate a, hence the value of the state variable 
Zoa in the microstate a, does not vary with time, the average Zo(t) for the macrostates 
varies with t: 

Zo(t) = EPa(t)Zoa. (22) 

a 

It should be pointed out that the entropy, as formulated in Eq. fl2T]) is also an average of 



In p) 



151 ]. the negative of the index of probability In p. There will be times, when we will 



also use an overbar such as in Zo(t) to indicate such averages for the sake of clarity. For 
common thermodynamic quantities such as average energy, volume, etc. the normal practice 
is to not use the overbar (unless clarity is needed) as it is mostly these average quantities 
that we deal with. 
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B. An Arbitrary Body 



It should also be stressed that the microstates for a body remain the same whether the 
body is isolated or not. We can apply Eqs. fl^T]) and (1221) to determine the entropy and 
the average quantity for any body, isolated or not (such as the system S or the medium E, 
neither of which is isolated). In the following, we will always use i to denote a microstate of 
a body but reserve a to denote the microstate of the isolated system. The entropy and the 
average energy of a macrostate of a body is given by 

S{t) = ~J2piit)\npi{t), (23a) 

i 

E{t) ^ EPiit)E,, (23b) 

i 

where i denotes one of its microstates, whose probability is denoted by Pi(t) > 0. While we 
can certainly allow microstates with probabilities Pi{t) = in Eq. (125]) . we find it convenient 
to only allow microstates with non-zero probabilities in the sum. Microstates with non-zero 



probabilities will be identified as allowed 



15| in this work. 



While there cannot be any doubt about the validity of Eq. (]23bp . one may feel some 



reservation about Eq. ( ]23ap for the entropy of an open system. Therefore, we will give a 
direct proof of Eq. fl23ap in Sect. El 

IV. INTERNAL EQUILIBRIUM THERMODYNAMICS 

A. Equiprobability Concept and Consequences: No Internal Variables 

We will first consider the case when there are no internal variables. As is the normal 
practice (see I for details), we assume that the medium is in internal equilibrium even if the 
system and the medium are not in equilibrium. This assumption is similar to the assumption 
of local equilibrium in the conventional nonequilibrium thermodynamics noted in Sect. Illli 
As discussed in I and in the review isj], the condition for the internal equilibrium to be met 
is that the entropy has the maximum possible value at each instant for the instantaneous 
average value Xie = ^(t) of the observable of the body. For the case considered in I, they 
are Eje = E{t) and Vie = V{t) for the medium; the number of particles of the medium is 
not allowed to change at all, so that A'ie = is always kept fixed. It is easy to see from 
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Gibbs' formulation of the entropy in Eqs. (12T]) or fl23ap that this happens if and only if all 
the microstates that are allowed {pi{t) > 0) at that instant are equiprobable. Another way 
to think about the internal equilibrium is to imagine isolating the medium by disconnecting 
it from the system. This will keep X fixed at Xie- Then the entropy of the isolated medium 
cannot change anymore. In other words, it is in equilibrium. We can apply the same idea 
to any body in internal equilibrium. The body will remain in equilibrium if it is isolated. 
Let us follow the consequences of this concept beyond what was discussed in I. 

(1) We allow for the possibility of external forces including stresses acting on the system 
and produced by the medium; see Fig. [1] Under internal equilibrium, 

p,{t) = l/Wit), Vz, (24) 
I — I 

where W{t) is the number of allowed microstates |15|] at that instant. This immediately 
leads to the Boltzmann entropy 

S{t) = lnW{t) (25) 
for a system in internal equilibrium, a common assurnption in non-equilibrium thermo- 



dynamics; see for example Bouchbinder and Langer {4 1 . Since the entropy is maximum 
at each instant t, it cannot increase further if all observables are kept fixed at their 
values Xie at that moment t, when the internal equilibrium is achieved. To ensure 
that X is held fixed at Xie, we isolate the system from the medium so that S turns 
into an isolated system. Let its entropy be denoted by S'is(t) after isolation. As the 
entropy is already at its maximum, it cannot change. In other words. Sis = 5'(Xie) at 
fixed Xie must be independent of time. 

The above argument proves that the entropy S'ie = S'(Xie) has no explicit t- 
dependence when the system is in internal equilibrium: 

^lE = 5(^iE, V^iE, ■■■ ,N) = S{E{t), V{t), ■■■ ,N). 

Its implicit time dependence when S is treated as an open system comes from the 
temporal evolution of X(t). Thus, the (maximum) entropy in Eq. (l25l) will change 
with time as X()f:) changes with time in the open system. For the isolated situation. 
Sis = 5'(Xie) at fixed Xie will remain constant. 
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Since the entropy is maximum for fixed Xje, there cannot be any additional irreversible 
entropy production diS^^^^ anymore 



diS'iElxiE = in internal equilibrium (IE); (26) 
here, the symbol l^j^ means that Xje are held fixed. 

When the system is in internal equilibrium, its various parts must be in equilibrium 
with each other. Otherwise, there would be irreversible entropy generation. 

It also follows from (3) that all the arguments that one uses to follow the consequences 
of equilibrium can be applied to different parts of S that are in equilibrium. For 
example, the arguments that establish that a system in equilibrium can only sustain 
uniform translation and rotation |l6, Sect. 10] can be applied without any change to 
a system in internal equilibrium. As this result is going to play an important role in 
our approach, we state it as a theorem. 

Theorem 1 There cannot be any relative motion between different parts ofT, for fixed 
Xffi in the state of internal equilibrium. Thus, a system in internal equilibrium can 



only sustain uniform translation and rotation ll d l. 



Proof. We refer the reader to Landau and Lifshitz [16l, Sect. 10] for the details. 
Their argument applies without any change to an isolated system in equilibrium. We 
easily extend their argument by considering our system at some instant t when it has 
XiE = X(t). We keep X(t) fixed at Xje by isolating the system from the medium. 
The arguments now apply without any change to the system in internal equilibrium. 
This proves the theorem. ■ 

The axis of the uniform rotation must be a principal axis of the instantaneous moment 
of inertia of the system. Otherwise, the system will undergo precession in space {26 1 
and the rotation will not be uniform. 

Even with internal equilibrium in the system, there are both elastic and inelastic or 
plastic deformations [271, Sect. 7.4.2], which result in viscoelasticity in the system. 



22 



If and only if the system is under internal equilibrium, the derivatives of S{t) with respect 
to X(t) have the significance of the field variables, also called the conjugate variables, which 
we will denote by y{t) or Y(t) : 

where X'(t) denotes all other elements of X(t) except the one used in the derivative; compare 
with Eq. ( 154|) . The temperature T(t) and the pressure P(t) are defined in the customary 
way by 

where X'(t) denotes all other elements of X(t) except E(t) or V{t), respectively in the two 
derivatives. The pair of quantities X(t), Y(t) or X(t), y{t) are called conjugate to each other. 

The definitions of the conjugate fields give us an alternative way to interpret the internal 
equilibrium. We imagine bringing the system in contact with another medium whose field 
variables are also Yje, where 

^ ^ (dsity\ 

Tie - \dMt)J X:.' 

here l^j^, denotes the value of the derivative at Xie- To distinguish this medium from 
the first medium that is characterized by Yq = {To,Pq,- ■■), we denote the first medium 
by S(Yo), and the new medium by S(Yie). The system S in internal equilibrium with 
observables Xje is in equilibrium with the new medium S(Yie). If the system is isolated 
by disconnecting it from S(Yie), then there cannot be any change in the macrostate of the 
system as all its observables remain constant at Xie. In other words, the system S after 
being isolated remains in equilibrium if it was originally in internal equilibrium, as noted 
earlier. This will not be true if there are internal variables, to which we will turn to in a 
moment. 



B. Zeroth and the Second Law 



All the above discussion has been for the entire system, but can be easily extended to 
a system consisting of various subsystems ak, each in internal equilibrium. Let us consider 
the system shown in Fig. [2|, which consists of two subsystems Ui and a2 whose internal 
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temperatures are Ti(t) and T2(t) > Ti(t), respectively. Let their respective energies be Ei(t) 
and E2{t), with their sum denoted by E{t). We consider all other observables fixed for both 
subsystems. We first consider the system to be isolated with no medium, as shown in Fig. 
Ulp). Then E remains constant, but not the individual energies. The irreversible entropy 
gain for the entire system is 



during an infinitesimal heat transfer dQ from the hotter subsystem to the colder subsystem. 
Here, 



plays the role of the thermodynamic force F^{t) for the isolated system, and should be 
contrasted with the same quantity for an open system in Eq. f|T8|) . As the system is 
isolated, this is also the entropy change S-^'^S = S^^Si + d^^^S2 for the system, with 

S^)s, = ^, S^^S2 



my ' T2it) 

for the two subsystems as follows from Eq. (fT5|) . We now consider the system to be in a 
medium at a fixed temperature Tq, as shown in Fig. |2t^a). We take Tq to be the equilibrium 
temperature of the isolated subsystems; it is intermediate between Ti(t) and T2{t). The 
infinitesimal heat given out by the hotter subsystem is now dQ = dQ' + dQ2. The heat 
gained dQ' + dQihj the colder subsystem must be exactly the heat loss dQ, since we are 
dealing with an isolated system. Therefore, dQi = dQ2, so that the entropy of the medium 
does not change. As the entropy change of the isolated system is equal to that of the system, 
we have 

since dQ = dQ' + dQi = dQ' + dQ2. This is the same irreversible entropy gain in Eq. 
fl29|) for the isolated system in Fig. [2](b). This should not be surprising as none of the 
heat transfers is isothermal. Thus, bringing the isolated system S, which consists of two 
subsystems, in contact with a medium, characterized by the equilibrium temperature Tq of 
S, does not affect the irreversible entropy production. It is easy to see that the arguments 
can be extended to many subsystems and to other field variables. We will not pause here to 
do that. 
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At this moment, it is important to follow another important consequence of the thermody- 
namic force -Fe(^), which vanishes if and only if the system has come to thermal equilibrium. 
This is the zeroth law of thermodynamics in terms of the internal temperatures of the two 
subsystems. Thus, the internal instantaneous temperature plays the role of a thermodynamic 
temperature in that the heat always flows from a hotter subsystem to a colder subsystem. 
The above result can be easily generalized to many subsystems. 

Let us now consider the volumes of the two subsystems to adjust as they come to equi- 
librium. All other extensive observables are considered fixed. Then the same reasoniong as 
above results in 

where dV{t) is change in the volume of ai; the volume of the isolated system remains 
unchanged. The corresponding thermodynamic force 

vanishes when the system comes to mechanical equilibrium. It usually happens that thermal 
equilibrium requires mechanical equilibrium in that the forces exerted on each other by any 
two subsystems must be equal and opposite. Thus, the conditions for the equilibrium is 
that not only the pressure P{t) but also the temperature T(t) have the same value in both 
subsystems. 

C. Presence of Internal Variables 

We can easily extend the above discussion to include internal variables by replacing 
X(t) by Z(t). In the context of internal variables, their conjugate variables are known as 
"affinity.''^ The general form of Eq. fl27p is 

W(t) fdS{t)\ 

where Z'(t) denotes all other elements of Z(t) except the one used in the derivative. The 
affinity a(t) is given by 

A(t) ^ (dS{t)\ 



so that w(t) consists of 
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and a(t). The generalization of the thermodynamic force in Eq. (1T8|) is given by 



Fz(t) = F [w{t)] = w{t) - Wo 



(34) 



with Wo representing the equihbrium value of w corresponding to the state variable Z. 

Let us now consider the system to be in internal equilibrium, while the medium containing 
it is S(Yo,Ao), where Yo,Ao characterize the medium. If we now disconnect the system 
from this medium, but bring it in contact with another medium S(Yie,Aie), where Yje 
and AiE are the field and affinity vectors associated with the system in internal equilibrium, 
then the system will remain in equilibrium with this medium. This is no different than what 
we have said above in the absence of any internal variable I. 

But the situation is very different when we try to keep the system isolated. Since internal 
variables are not under the control of the observer, they cannot be manipulated to remain 
constant after isolation and will continue to change. Thus, after the isolation, Xis, the value 
of X(t) at the instance of isolation, will remain constant, but l(t) will not remain fixed at 
its value lis at the instant it was isolated. This time variation in the internal variables is 
what drives this isolated system towards its equilibrium state during which its entropy will 
continuously increase. This is very different from the case above when there were no internal 
variables. Thus, a system in internal equilibrium cannot be isolated as an equilibrium system 
if there are internal variables present. It can only remain in equilibrium with the medium 



D. Maxwell Relations 

As the concept of internal equilibrium is not that different from the concept of equilibrium, 
it should not come as a surprise that there are analogs of Maxwell relations. We recall that 
in equilibrium thermodynamics, the standard Maxwell relations for a system characterized 
by only S and V (fixed A^) QbTQ clS follows in terms of Jacobians: 



All four Maxwell relations use the same numerators 9(To, S, N) and d{Po, V, N). They use 
different denominators. Thus, they can all be combined into one compact relation that can 



S(Yis,Ais). 



d{To,S,N) _ d{Po,V,N) d{To,S,N) _ d{Po,V,N) 

diV, S, N) diV, S, N) ' diPo, S, N) d{Po, S, N) ' 

d{Po,V,N) ^ din,S,N) d{To,S,N) ^ diPo,V,N) 

9(To,y,iV) 9(To,y,iV)' diPo,To,N) diPo,To,N)- 



(35) 



26 



be simply written as 

d{To,S,N) = d{Po,V,N). (36) 

Here, the relation only has a meaning if each side is divided by one of the possible denomi- 
nators d{V, S, N), d{Po, S, N), (9(To, V, A^) and d{Po, Tq, A^) on both sides. 

We now consider a system. For simplicity, we assume that only one internal variable, 
which we denote by ^, characterizes this system. We assume the system is in internal 
equilibrium. To simplify the notation, we will suppress but use the additional variable 
^ along with the other two variables. By considering the system in a medium S(Yis, Ajs), 
we recognize that the system is in equilibrium. Thus, the Maxwell relations now for fixed ^ 
can be compactly represented by 

diT,S,0 = diP,V,0, (37) 

by replacing Tq, Pq by Tit) = Tjs, P{t) = Pis, which for simplicity have been written as T, P. 
The extension to the Maxwell relation in terms of the internal variable requires considering 
the pair A, ^ in place of T, S or P, V, where A denotes the conjugate affinity to ^. For fixed 
V, the Maxwell relation is 

d{T,S,V)=d{A,tV), (38) 
and for fixed S, the Maxwell relation is 

d{P,V,S) = -d{A^,S). (39) 

If we consider the system in the medium E(Yis, Ais), then we have the standard equilibrium 
Maxwell relations similar to those in Eq. (1351) . However, we are interested in the possible 
"non-equilibrium" Maxwell relations when the system is in the medium S(Yo, Aq). We first 
consider fixed C,- The Maxwell relation in Eq. fl37|) turns into the identity 



d{T,S,0 _ d{P,V,0 
d{Po,S,0 diPo,S,0 

that is 

dT\ fdP\ fdV\ _(^\ 



for the enthalpy and 

d{T,S,0 _ d{P,V,0 
d{To,V,0 d{To,V,0' 
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(40) 



(41) 



that is 



for the Helmholtz free energy; the details will be given in a separate publications [28[. One 
can also obtain Maxwell relations at fixed V or 5*. For example, we find the following 
Maxwell relation 

d{P,V,S) _ d{A,^,S) 
d{^,V,S) d{^,V,S)' 

that is 

fdP\ _ fdA\ 



Similarly, from the Maxwell relation 

diT,S,V) _ d{A,^,V) 
di^,S,V) di^,S,V)' 

we find 

fdT\ _ _ rdA\ 
\'d^)v,s~~\ds),y- 
The Maxwell relations in Eq. ( I371l39l) contain the internal fields of the system and not of 



the medium 



when the system is out of equilibrium. Obviously, the extensive variables 



in the relation must refer to the system. 



E. Internal Equilibrium Thermodynamics versus Local Thermodynamics 

We will argue in Sect. |V]that the concept of internal equilibrium, which we adopt, is no 
different than the concept of local equilibrium used in Eq. (|5]) or in the Gibbs fundamental 
relation in Eq. ([7]). Despite this, the two approaches based on the concept of local and 
internal equilibrium, respectively, differ in some important ways that will be elaborated 
later. Here, we briefly mention these differences. The first important difference is that our 
approach is truly a statistical mechanical approach for non-equilibrium systems. Once the 
probabilities of microstates are known, the averages and other moments of all state variables 
that are used to identify the microstates and the entropy are determined for the macrostate. 
For example, the average fluctuation in Z for a body is given by 

[AZ{t)f^^p,{t)[Z,-Zit)]\ (42) 
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where 

m = Ep.WZ. (43) 

i 

is the average Z for the body; compare with Eq. (122|) . As the microstate probabihties exist 
even when the system is out of equihbrium, these averages including the entropy exist at 
all times even if the system is not in internal equilibrium. Their temporal variations are 
controlled by the dynamics governing the system and give rise to various balance equations; 
see for example Eq. f fT2|) . The use of a probabilistic approach in the determination of 
the entropy and other statistical properties means that the dynamics in the system must be 
stochastic and not deterministic, as the entropy remains constant in a deterministic dynamics 
15| . A consequence of the stochastic nature is that irreversible dissipation becomes an 
integral part of the statistical description of any system, which then results in the law of 
increase of entropy as captured by Eq. ([2]). It is this law that was the foundation of our 
approach in I, and which we develop further in this work. 

The second difference is in the identification of the thermodynamic potentials for the 
open system, which has been discussed in I and will be further elaborated in Sect. IVII and 
again in Sect. IX B[ It is discovered in our approach that thermodynamic potentials con- 
tain the field parameters (temperature, pressure, chemical potentials,etc.) of the medium, 
which determine how far an open system is from its equilibrium with the medium. These 
thermodynamic potentials satisfy the second law in that they do not increase in a sponta- 
neous process. In contrast, the form given in Eq. (j6]) for the local Gibbs free energy density 
or its integral over the volume of the system does not always satisfy this requirement; see 
the discussion surrounding G{t) in Eq. (22) of I. However, there is no discrepancy for the 
internal energy in the two approaches as both approaches give the same Gibbs fundamental 
relation. This is because the fundamental relation in both approaches does not depend on 
the field parameters of the medium, but include the instantaneous field parameters of the 
system. 

The third difference is in the reversible and irreversible entropy changes, which depend 
on the equilibrium value wq {yo or qq) of the conjugate field w {y or a), as is easily seen from 
Eqs. (|T6l) . ( fT71) . (|T9l) and (1201) . In the localthermodynamics, these quantities are determined 
by the local conjugate fields. 

The presence of the medium field variable in the thermodynamic potentials of the system 
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does not mean that the situation would be any different if we consider the system or sub- 
systems to form the isolated system Sq without any medium, as shown in Fig. Et^b). In this 
case, which we consider in Sect. IXIBt we again find that the thermodynamic potentials do 
not depend on the field variables of the system or subsystems. The role of the field variables 
of the medium S(Yo, Aq) are now played by the equilibrium conjugate variables Yq, Aq of 
the system or subsystems. This then leads us to the following important theorem: 

Theorem 2 An isolated system S is no different than the open system S in an extensively 
large medium S(Yo, Aq), provided the medium is appropriately chosen to represent the equi- 
librium state (in terms of Yo,Aq) of the isolated system S. 

In particular, the reversible entropy change and the irreversible entropy generation in the 
two cases (a) and (b) in Fig. [2] are exactly the same, as they both depend on the equilibrium 
conjugate variables Yq, Aq of the system or subsystems. An example of this is already seen 
in Eqs. (fT6ll and (ITTl) . We defer the proof of this theorem to Sect. IXIBI 



V. additivity of entropy and QUASI-INDEPENDENCE 

For simplicity of discussion, we consider all systems to be stationary in this section, so 
that we only deal with internal energies. It was noted in Sect. IIIII and recently reviewed 
in the entropy of an isolated body is given by the Gibbs formulation in Eq. fl?Il) . 
regardless of whether it is in equilibrium or not. There is no reason to believe that this 
formulation also applies to an open body under all conditions, though its applicability in 



equilibrium is not in dispute 16|]. We now prove that this formulation also applies to an 
open body under a condition that is always taken for granted. We will specifically consider 
our system S at some instant t, but the conclusion is valid for all bodies. Let us consider 
all allowed microstates of S with fixed number of particles N; we index these microstates 
by i = 1, 2, ■ ■ ■ ,W{t). Theses microstates correspond to all possible energies and volumes 
of the system. We use 5 to denote the microstates of S whose number of particles N is 
also fixed. A specification of the microstates i and 5 gives a unique microstate specification 
a representing a microstate of the isolated system Eq. Hence, the number of microstates 
Wo{t) of the So is the product 

Wo{t) = W{t)W{t), (44) 



30 



where W{t) and W{t) are respectively the number of all allowed microstates 15|] of E and 
S, respectively at that instant t. 

As the concept of microstates does not depend on the nature of interactions (they exist 
even in the absence of interaction), the above equation is valid for all kinds of interactions. 
Let Eo, E{t) and E{t) denote the internal energies of Sq, S and S, respectively. Let i?o^'^*^(t) 
denote the mutual interaction energy between S and S at that instant. For short-ranged 
interactions, this energy is determined by the surface dV{t) of S. For convenience, we 
assume that this entire area is exposed to the surrounding medium, as shown in Fig. [H If 
long-ranged interactions are also present, or if the system size is very small, this energy may 
depend on the entire volume V{t) of S. In all cases, this energy is defined by the following 
identity 

E^ = E{t) + E{t) + E'i'''\t). (45) 
Because of the smallness of S relative to Sq, E{t) « E{t). If it happens that 

Et'\t)«E{tl (46) 

we call S and S quasi-independent. For quasi-independence, the linear size of the system 
must be at least as large as, but hopefully larger than, the correlation length in the system. 
In this case, we can neglect their mutual interactions, which is a common practice in the 



discipline [l6|. The quasi-independence of the system and the medium holds to a very 
high degree of accuracy for all short-ranged interactions provided the system itself is 
macroscopically large so that the ratio of its surface to volume is insignificant. In most 
cases, this will also ensure that the correlation length is small compared to the size of the 
system. If there are also long-ranged interactions, then we can still have quasi-independence 
provided these interactions are relatively weak and shielding occurs and that Eq. (l46l) and 
the condition on the correlation length hold simultaneously. 

Let us now assume that S and S quasi- independent. In this case, the microstates of the 
two systems are independent of each other to a very high degree of accuracy and we have 
(we suppress all state variables for simplicity of notation) 

Pa{t) = Pi{t)Pa{t)- 

Now, using 

hipi{t)pa{t) = Inpi(t) + \nps{t), 
31 



and the sum rule 

i:p5(t) = i, ep.w = i, 

a i 

we find that 

i a 

where the two terms in the above equations represent the entropies of the system and the 
medium 

S{t) = -EP^it) Inp. Sit) = -ZPait) Inps(t), (47) 

i a 

respectively. This demonstration justifies the additivity of entropies 

Soit) = Sit) + Sit) (48) 

as a consequence of quasi-independence so that 

Eo = Eit) + Eit) (49) 

also holds to a very good approximation. Note that we have neither assumed the medium 
nor the system to be in internal equilibrium in the above demonstration. 

If the system and the medium fail to be quasi-independent because their mutual inter- 
action cannot be neglected, then Eq. fH^ is violated. In this case, the presence of this 
interaction acts as a constraint on Sq. Consequently, the entropy now will be strictly less 
than the above entropy in Eq. fHHl) . We denote this difference by 5'o"*^(t) < 2^, which is 
defined by the following identity 

S,it)^Sit) + Sit) + st'\t). (50) 

This identity reduces to the above additivity in Eq. ( HHl) provided 

St'\t) «Sit); (51) 



compare with Eq. fl46l) . This inequality will in general hold only if the interaction energy is 
also negligible. 

If the strong inequality in Eq. ( l5Ti) is not satisfied, we have lost the additivity property 
of the entropy. Let us assume that the strong inequality is satisfied for some large size of 
the system S. Now, as the size of the system decreases, which is what will happen on the 
way to considering physically infinitesimal volume elements used in Eq. ([3]), there comes a 
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point where the strong inequahties in Eqs. (H6|l and (15T|) are violated. This will destroy the 
additivity of the entropy as exhibited in Eq. f l48|) . 

The violation of entropy additivity occurs at intermediate sizes of the system, somewhere 
between the macroscopic size where Eqs. fHB]) and fl^T]) are valid, and small local or mi- 
croscopic size containing a small number dN of particles. For example, for dN ^ 10^*, 
the surface to volume ratio for the volume element dV is about 10~^, implying an almost 
imperceptible error in neglecting the interaction entropy SQ^^\t), provided the linear size 
of this region is large compared to not only to the inter-particle separation [2, p.l] but 
also the correlation length in the system. Under these conditions, the integrand in Eq. 
truly refers to a "physically" infinitesimal volume element containing a very large number of 
particles. In this sense, our startirig premise is similar to that adopted in the conventional 
non-equilibrium thermodynamics |9l-|l2|. In particular, our concept of internal equilibrium 
is no different than the concept of local equilibrium in conventional non-equilibrium thermo- 
dynamics, as observed in I, except that we require quasi-independence, which imposes the 
strong condition that not only the interaction energy be small but also that the linear size 
be larger than the correlation length. At present, there is some evidence that the correlation 



length in a g 



transition 



30|. 



ass forming system appears to increasing as the system approaches the glass 



The discussion above also clarifies that the additivity of entropy is a consequence of the 
additivity of energy of various parts of the body and that the interaction energies between 
them must be negligible. This additivity of the energy and entropy was adopted in I. 
However, as we will be interested in considering parts of S as subsystems in this work, the 
additivity of their entropy requires that their mutual interaction energies be also negligibly 
small compared to their individual internal energies, and that their linear sizes be large 
compared to the correlation lengths 30|. These requirements put a strong condition on the 
sizes of subsystems. 
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VI. THERMODYNAMIC POTENTIALS 



A. Fixed Number of Particles N of the System S 

Despite similarities between our approach and that adopted in conventional thermody- 
namics IJQ, there were important differences noted in I. One of these was the discovery 
that the differences of the temperature and pressure of S and S, which are the same as that 
of S and Sq, play the role similar to internal variables. The second difference was that the 
Gibbs free energy in our approach exists even if the system is not in internal equilibrium, and 
involves the temperature and pressure of So or S; of course, we assume that the medium is in 
internal equilibrium; see I and Sect. ??. Under this very weak assumption for the medium, 
its field variables such as the temperature To, pressure Pq, etc. are well defined, and are 
unaffected by whatever processes happen to be going on within the system S or whether E 
is homogeneous or inhomogeneous. When the number of particles N in the system S is held 
fixed, the appropriate thermodynamic potential is the Gibbs free energy, which is identified 
as 

G{To,Po,t) = E{t)-ToS{t) + PoV{t), (52) 

where the observables E{t),S{t) and V{t) have explicit time-dependence for fixed N. The 
particular form of the Gibbs free energy is in accordance with the second law in Eq. (j2]), 
and remains valid even if the system is so far out of equilibrium that its temperature and 
pressure cannot be defined; see also Landau and Lifshitz see Sect. 20]. It is not 
surprising, therefore, that it contains the temperature and pressure of the medium, which 
are well defined. Theses quantities do not exists for the system unless it happens to be at 
least in internal equilibrium. Even then, the Gibbs free energy is given by Eq. fl52l) and 
contains the temperature and pressure of S. As such, it does not represent a state function 
of the system. In engineering context, this quantity is also known as exergy or availability 



m 



However, our form in Eq. ( I52l) differs from the local form of the Gibbs free energy in Eq. 
dH]), which contains the local temperature and pressure. One can argue that the identification 
of the Gibbs free energy in I was for the entire system, but that once we account for the 
inhomogeneity by considering subsystems, the Gibbs free energy for each subsystem will 
somehow become consistent with that in Eq. (Q. This is a reasonable possibility and we 
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need to investigate this possibility. This issue is deferred to Sect. IX B[ 

As the method to identify the Gibbs free energy or other thermodynamic potentials that 
follow from the second law in Eq. is going to be employed here several times, we briefly 
sketch the derivation for the sake of continuity. The full details are given in I. We do not 
assume the existence of the temperature, pressure, etc. of the system to include the situation 
in our discussion when the system is far away form equilibrium so that they are not defined. 
For simplicity, we consider a monatomic system of structureless particles under no external 



shear as in the previous work [1]. Accordingly, we only consider the energy E, volume V 
and the number of particles to describe the macrostate of the system at any instant t. 
No internal variables will be considered at this moment. The system and the medium are 
assumed not to be in equilibrium. 

We use the additivity in Eq. fHHj) to write the entropy So(t) of Sq as the sum of the 
entropies S{t) of the system and S{t) of the medium. If we also assume that the latter is in 
internal equilibrium, then we have 



there is no explicit t-dependence in S{E, V, N) due to its internal equilibrium. With A^ and 
A^ fixed, we expand 5*0 in terms of the small quantities E{t) and V{t) of the system 



where the derivatives are evaluated at Eo,Vo,N. However, as A^ is very close to A'o, there 
is no harm in evaluating the derivatives at Eq,Vq,Nq. This is the reason that we have used 
Xo above for the derivative. This approximation will be made throughout in this work. 
The error is inconsequential when the system is a very small part of the isolated system. It 
follows from the internal equilibrium of S that 



We observe that S = S{Eq, Vq, N) is a constant, which is independent of the system E. 
Thus, 



n 



So{Eo, Vo, No, t) = SiE, V, N, t) + S{E, V, N); 



(53) 





(54) 



So{t) - 5 ~ S{E, V, N, t) - E{t)/To - PoV{t)/To. 



(55) 
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In terms of 

G{t) = H{t) - ToS{t), H{t) = E{t) + PoV{t), (56) 

we finally have 

Soit) -S = Sit) - Hit) /To = -Git) /To. (57) 

It immediately follows from Eq. ([2]) that the Gibbs free energy Git) of the system in Eq. 
( I52l) continuously decreases as the system relaxes towards equilibrium, a result quite well 



known in classical thermodynamics 16|: 



'-^ < 0. (58) 

dt - ^ ^ 

The function Git) continues to decrease and finally becomes identical to the equilibrium 
Gibbs free energy at the current temperature and pressure To, Pq- If we abruptly change the 
temperature T' and pressure P' of the system in some state A', where the system was in 
equilibrium, to a new state A where the temperature and pressure are Tq, Pq? respectively, at 
time t = 0, then the initial values of the energy, volume and entropy at the new temperature 
and pressure remain equal to their respective equilibrium values in the previous state A' as 
the microstate probabilities Pait) at t = have not had any time to change. Thus, initially 

G(0) = E'^- ToS'^ + PoVi 

in the state A; the quantities with a prime are the equilibrium values in the state A'. The 
Gibbs free energy decreases in accordance with Eq. (!58l) and eventually becomes equal to 
its new equilibrium value 

Ga = Ea — ToSa + PqVa, 

where the quantities with the subscript A denote the equilibrium values in the new state A. 

It should be noted that the equilibrium Gibbs free energy in the state A' before the 
abrupt change is 

G'^ = E'^- rs'^ + P'Vi, 

so that the Gibbs free energy undergoes a discontinuity at t = due to the abrupt change: 

AG'^ = iT'-To)S'^-iP'-Po)Vi. 

Its magnitude and sign has nothing to do with the second law as the abrupt change is not 
a spontaneous process. 
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A similar looking quantity G(t), see Eq. ([6]) for its local analog in the local non- 
equilibrium thermodynamics |9l-|l2|. 



G{t) = H{t) - T{t)S{t), H{t) = E{t) + P{t)V{t), 



(59) 



which can be defined only when t 
was shown to increase with time 



le system is under internal equilibrium and not otherwise, 
l| during relaxation 



dG{t) 
dt 



> 



(60) 



in a cooling process. Since it does not always decrease with time, it cannot be taken as 
the Gibbs free energy; the latter is supposed to never increase as the system equilibrates 
spontaneously as happens with G{t); see Eq. (l58l) . 



B. Fixed Volume V of the System S 

Instead of keeping the number of particles in S fixed, let us keep its volume V fixed so 
that the volume of the medium is also kept fixed. The number of particles of the medium 
is no longer fixed. The entropy S{E, V, N) of the medium in Eq. fl53|) is expanded in terms 
of small quantities E and of the system. We follow the steps similar to those above and 
obtain _ _ 

~ ~ ~ / r) Q \ / r) Q \ 

N{t), (61) 



SiE,V,N)^SiEo,V,No)- | ^ 



E{t) - 

^ dN 



where Xq stands for Eq^Vq^Nq for reasons explained above in deriving Eq. (IMj) . Let us now 
introduce the chemical potential yUo of the particle in the medium by the standard definition 



dS 
dN 



Xo 



To' 



We thus find that in terms of 5 = S{Eq, V, Nq) 

So{t) - 5 ~ S{E, V, N, t) - E{t)/To + fioN{t)/To = -[E{t) - ToS{t) - fioN{t)]/To, (62) 
which identifies a different thermodynamic potential in this case as 



nit) = Eit)-nSit)-^oNit); 



(63) 
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this thermodynamic potential also uses the field variables of the medium. We should empha- 
size again that no assumption about the internal equilibrium of the system has been made. 
The system may or may not be in internal equilibrium. The application of the second law 
in Eq. ([2]) now gives 

^ < 0. (64) 

dt - ^ ^ 

C. Fixed N and V of the System S 

If both and V are kept fixed, it is easy to follow the above derivation to conclude that 
the Helmholtz free energy 

F{t) = E{t) - ToS{t) 

must continuously decrease as the system reaches equilibrium: 

dF{t) 



dt 



< 0. 



Thus, the second law for an open system is expressed in terms of different thermodynamic 
potentials depending on which variables are held fixed. 



VII. INTERNAL VARIABLES 



As said above, a suitable equilibrium macrostate description of the system requires a set 
of independent macroscopic observables that can be controlled by an experimentalist and 
whose values will allow the experimentalist to differentiate between different macrostates 
of the same system. It normally happens that experimentalists have a far less number 
of external controls than the possible extensive variables that can be used to characterize 
the macrostates. Thus, one does not characterize a macrostate, especially an equilibrium 
macrostate, by specifying all of the relevant extensive system quantities. For example, for a 
single component system, one normally uses E, V and to specify the macrostate if there 
are no external shearing forces. Let us for the moment consider a system without exter- 
nal shear. Usually, one considers a system with fixed A^; then E and V can be controlled 
by the two external variables Tq and Pq associated with the medium. However, these ex- 
ternal variables need not necessarily control the local or internal structures in the system 
at all times during its evolution towards equilibrium. As Frenkel has observed, the local 
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structures can be important when considering the structural relaxation in a glass or other 



non-equilibrium systems 
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p. 208]. For example, one can consider the average numbers 
of neighbors and next-neighbors of a given particle to describe the local structure in the 
system. These quantities multiplied by the number can play the role of internal vari- 
ables. The corresponding conjugate variables, normally identified as "chemical potentials" 
or "affinity" for these internal variables usually vanish in equilibrium. Frenkel goes on and 
calculates viscoelastic effects due to structural changes and compares them with Maxwell's 
model of elastic relaxation or an RC-circuit. This investigation by Frenkel [3^ shows that 
internal variables can play an important role in the temporal evolution towards equilibrium 
in some systems such as glasses. As such, they become an integral part of the description of 
any non-equilibrium system and determine the relaxation of the system Sect. 78]. The 
internal variables are also called hidden variables or internal order parameters. 

To introduce the concept of internal variables, let us consider our isolated system Eg 
for which one can identify a set of conserved quantities, i.e. integrals of motion. For a 
mechanical system of s degrees of freedom, the number of such integrals of motion are 
2s — 1 261]. Of these integrals of motion, those that are additive play an important role in 
thermodynamics and statistical mechanics. The notable ones are the energy, and linear and 
angular momenta of the system, among others such as the polarization, magnetization, etc. 
For the moment, let us consider Eg to be stationary. Its macrostate A^o is characterized by 
fixed internal energy Eq, volume Vq, particle number A^o and other extensive observables, 
collectively denoted by Xq. All these observables are constant for Sq. Let us consider the 
energy Eq, which is an integral of motion. It usually happens (see below for an example) 

(k) 

that there are many different components Eq of the energy whose total sum is the energy 
of Eq: 

n+l 

Eo^J2^o\t), (65) 

k=l 

where n + l > 1 is the number of energy components. It is Eq that is a constant of 
motion, not the individual components E^'^^t); the latter will continue to change as the 
system evolves in time while maintaining Eq. (1651) . Let Wo(Ko,t) denote the number of 
microstates corresponding to the macrostate at time t. At each instant t, the microstates 
in Wo(X.o,t) can be partitioned into groups according to the possible values of EQ^\t). 
Because of the sum rule in Eq. ( |65l) . only n of the components are independent for a given 
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Eq, which we take to be given by k = 1, 2, ■ ■ ■ ,n. We will denote this set by an n- vector 
Io(t) whose elements are| Jq'^'* = EQ''\t), k = 1,2, ■ ■ ■ , n|. Then, 

n 
k=l 

Let Wo(KQ,lQ{t),t) denote the number of microstates for a given Zo(t). These microstates 
define a new macrostate, which we denote by Aq. Obviously, 

Wo{Xo,t) = 5^W^o(Xo,Io(t),t); (66) 

the sum is over all possible Io(i)- As the system evolves, different components i?o'^^(t) of lo(^) 
evolve in time t, but Eq and Xq remain fixed. Thus, a better understanding of the evolution 
of the system can be obtained by monitoring how the various components E^''^ (t) change in 
time. For this, it is better to use Zq = (Xq, Io(^)) to identify the macrostate Mq even though 
individual E^\t) cannot be controlled by the observer. As Eq can be controlled by the 
observer, it is still the choice observable to be used for identifying a macrostate. This is even 
more true for the isolated system for which Eq is a constant of motion. The n components of 
Io(t) then play the role of internal variables in developing non-equilibrium thermodynamics 
of the isolated system. 

It is evident that the same "extended" description should also be useful for an open 
system S. The only difference between an open system and an isolated system is that not 
all elements of X remain fixed. Some of the observables, denoted by X' are controlled by 
external field parameters Yq (such as Tq,Pq, etc.) of the medium so that they do not remain 
fixed but continue to fluctuate about their mean X'(t) that keeps changing in time. However, 
at least one of the extensive observables such as N must be kept constant to quantify the 
size of the system |33|. Thus, for an open system, these observables can be replaced by 
the fields Yq, with the remaining observables remaining constant. We will denote the latter 
observables by C to remind us that they are constant. The open system can be either 
specified by X'(t),C or Yq, C. However, for the sake of convenience, we will continue to 
use X(t) rather than X'(t),C or Yq, C. Let us now consider S which is not in internal 
equilibrium so that it undergoes internal deformation due to relative motions between its 
various parts. If there are external strains on the system, they can be controlled by us from 
the outside. Hence, they will not be considered as internal variables. However, internal 
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stresses acting on various parts of the system when there are no external strains on the 
system are beyond our control and must be treated as internal variables in describing the 
system. As said earlier, we can describe the internal forces acting on each part in terms 
of translation and rotation of its various parts; see Sect. Ill A[ These motions must be 
described by the use of suitable internal variables such as the linear and angular momenta, 
as was discussed in Sect. II B[ 

As the internal variables are uncontrollable, their affinity in equilibrium must vanish as 
we prove now as a theorem. 

Theorem 3 The affinity of an internal variable must vanish in equilibrium. 

Proof. It is sufficient to prove the theorem for an isolated system. Also, we will prove it for 
the energy components in Eq. (jHSD- The extension to the general case is a trivial extension 
and will not be done here. As we are dealing with equilibrium, we consider equilibrium 
values of all the quantities, which are going to be represented by suppressing the argument 
t as they are stationary. We now construct the following partition function for the isolated 
system 



where ao is the the n- vector |ao | equilibrium affinity. Such a partition function correctly 
describes the situation in which the n internal variables are not constant but keep changing 
from microstate to microstate. We now observe that this partition function reduces to the 
equilibrium value Wq(X.o) in Eq. fl^ (where we take the limit t — )■ oo), provided A^^'^ = 0: 



for each of the internal variable in the set Io,eq- This proves the theorem. ■ 

The above theorem deals with equilibrium affinities, and says nothing about the affinities 

of the internal variables when the system is out of equilibrium. 

Let a denote one of the microstates associated with the macrostate Aio, and /3 one of the 

microstates associated with the macrostate Ao- Then, using their probabilities Pa(t) and 

P/3{t), we can determine the entropy of the two macrostates using the Gibbs formulation in 





-Aq ^ 0, k — 1, 2, ■ • • 



n, 
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Eq. m- 



'S'o(Xo, t) 



(67a) 



a 



/3 



(67b) 



For a macroscopically large system, the following standard statistical mechanical arguments 
can be used to highlight the maximum of the summand in Eq. (I66p . Let the maximum of 
the summand be denoted by Mo(t), which occurs for some particular value To(t) of Io(t) : 



We separate the maximum contribution from the sum and rewrite Eq. fl^B]) as follows: 



where the sum is over all remaining Io{t). It is normally the case that the ratio in the above 
sum is vanishingly small for a macroscopic system and that the sum can be neglected. In 
this case, we have 



For a macroscopically large open system such as S, the above equation is formally valid, 
except that we must replace Xq by X(t), which stands for X'(t), C, and Io(t) by I(t): 



We now prove an important theorem about the nature of the entropy. 

Theorem 4 The entropy expressed only in terms of the observables when (independent) 
internal variables are present must explicitly depend on t. 

Proof. We first consider the isolated system Sq and prove the theorem for it. For lo{t) 
to be independent of (fixed) Xq, it must surely have an explicit dependence on time. In 
other words, lo{t) must be a function of Xq and t. Let us assume that there is an explicit 
i-dependence in both Wo-functions in Eq. ( |68l) . As the entropy of the macrostate A^o is 
given by the sum over all microstates Wo(X.Q,t) in Eq. fl23ap . it must explicitly depend on 
t. Thus, the theorem is satisfied. If, however, neither of the VFo-functions in Eq. (|68|1 have 



Mo{t) = Wo{Xo,Mt),t) 




Wo{Xo,t)^Wo{XoMt),t). 



(68) 



W{X{t),t)^W{X{t),I{t),t). 



(69) 
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any explicit t-dependence, then this is possible only if To{t) becomes a function of Xq as the 
left side is only a function of Xq. Since Xq is constant, lo(t) itself must be constant. The 
latter is the situation in equilibrium: 

Io,eq = lo(Xo) a constant. (70) 

It follows that Iq = Io,eq is no longer an independent variable when the system is in 
equilibrium. Obviously, this case is not covered by the theorem since Iq is not independent. 
The entropy in this case is given by the Boltzmann formulation, cf. Eq. ( !25l) . and we have 
fromEq. 

S'o(Xo) ~ S'o(Xo,Io) a constant. 

Let us now consider the special case when the macrostate Aq satisfies the condition of 
internal equilibrium. In this case, WqCKq, Io(t), t) does not explicitly depend on t and should 
be written as Wo(X.Q,lQ(t)) with To(t) having an explicit time dependence. This entropy is 
again given by Eq. ( !25l) : 

5o(Xo,Io(t)) = InWoiXoMt))- (71) 

It now follows from Eq. ( pSj) that Wo(X.Q,t) must have an explicit time-dependence due to 
the explicit t-dependence of Io(t) in Wo(X.Q,lQ(t)). This is because different values of Io(t) 
will result in different values of PVo(Xo, Io(^)), which can be treated as Wo(X.o,t) associated 
with the macrostate at different times. This is the first case considered above. Thus, 
S'o(Xo,t) will have an exphcit t-dependence even though S'o(Xo, Io(t)) does not. 
This proves the theorem for an isolated system. 

Let us consider an open system such as S. Again, T(t) must be a function of X()f:) and t 
to remain independent of X(t). Let us assume that there is an explicit t-dependence in the 
ly-functions in Eq. (169|) . As the entropy of the macrostate of S is given by the sum over 
all microstates W(K{t),t) in Eq. fl67ap . it must explicitly depend on t. If, however, neither 
of the ly-functions have any explicit t-dependence, then I{t) becomes a function of X(t). In 
this case, it is not independent of X(t). This situation is then not relevant for the theorem. 

The possibility in which I(t) is independent of X(t), but W (X.{t) ,T{t)) has no exphcit 
t-dependence, when the system is under internal equilibrium, is very important. Fixing 
Xig = X(t) allows us to think of the system as an isolated system. Now, we can use the 
argument given above for the isolated system to conclude that different values of I(t) will 
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result in different values of iy(Xis, which can be treated as W{X.is,t) associated with 
the macrostate A4 at different times. In other words, the macrostate Ai does not represent 
an internal equilibrium state. Thus, we conclude that a macrostate J\f under internal 
equilibrium results in a macrostate Ai; the latter is, however, not in internal equilibrium. 
This proves the theorem. ■ 

It follows from the above discussion that a general thermodynamic state can be taken to 
be a function of internal variables along with other observables and time t when we deal with 
non-equilibrium states. For an open system in which many of the observables are controlled 
by external field parameters Yq (such as Tq,Po, etc.) of the medium, we can express leq 
either as 

leq = I(Xgq, C), 

or as 

Ieq^T(Y[„C). 

Away from equilibrium, the internal variable I(X'(t), C) differs from its equilibrium values 
Igq, and is normally treated as an independent variable and plays an important role in the 
dynamics of the system as the latter strives to reach equilibrium. Thus, it is not surprising 
that internal variables are employed to specify the macrostate of a glass. In non-equilibrium 
thermodynamics, this fact has been recognized for quite some time jol I2I1. 

Internal variables can also be related to the presence of internal degrees of freedom in the 
particles of interest. The internal degrees are more common in polymers but can also occur in 
small molecules in the form of rotation about some internal axes. An example will clarify the 
point much better. Consider a polymerization process resulting in a system of polydisperse 
linear polymer chains of average molecular weight M in a solution 3^. The model is defined 
on a lattice of sites and volume V = Nvq, with vq a constant representing the volume 
occupied by a lattice site. One normally uses E, V, M defined below in Eq. fl72l) . and the 
number of chains p as the standard observables that can be used to identify the macrostate 
(equilibrium or not) of the polymer solution. In turn, these quantities are controlled by the 
temperature, pressure and the initiation-termination and propagation rates; the last two can 
be related to the initiation-termination activity controlling the number of endgroups, two for 
each polymer, and the middlegroup activity. These activities determine the corresponding 
affinity or "chemical potentials." Let A^m = N — denote the number of monomers, each 
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monomer occupying a lattice site, in terms of the number of voids or sites not covered by 
monomers A^v so that 

TV 

M=— . (72) 

V 

In terms of the number of middle groups A^m = A^m — 2p, or A^m, the number of chemical 
bonds in the p polymers is given by 



There are two kinds of energy in the model 3J]. One kind of energy is due to mutual 
interactions of voids (v) with the end (E) and middle (M) groups, and the mutual interac- 
tions between chemically unbonded M and E. Let Nij,i,j =v,M or E, denote the number of 
nearest-neighbor contacts ij,i 7^ j, and Eij the corresponding interaction energies, respec- 
tively. The other kind of energy is due to intrachain gauche bonds (g), and hairpin turns 
(hp). Their energies are Eg = EgNg for gauche bonds and i?hp = f^hpA^hp for hairpin turns; 
here A"g, and A^hp denote the number of gauche bonds and hairpin turns and parallel bonds 
and eg, and e^p are their energies. In addition, there is a mutual interaction energy between 
two parallel (chemical) bonds, which may belong to the same or different polymers. Let 
A^P denote the number parallel bonds, each of energy ep. Thus, 

E= J2 £ijNij + egNg + ekpNkp + epNp= ^ Eij + Eg + E^p + Ep, (73) 

where we have introduced Eij, Eg, Ei,p and Ep with obvious definitions. We thus observe 
that the energy can be partitioned into six extensive energies, five of which can be taken as 
internal variables. 

To summarize, we conclude that the quantities that cannot be controlled by the observer 
can be identified as the internal variables. This statement should not be taken literary as 
what is considered uncontrollable today may not remain so in the future. Thus, to some 
degree, the decision to identify the internal variables is left to the observer. For us, any 
variable that cannot be controlled to have a fixed value when the system is out of equilibrium 



will be taken as an internal variable 35|. It should also be noted that the number of internal 
variables is not a unique number for a given system. For example, to describe local structures 
in a monatomic system [32] , one can consider any number of neighboring particles (neighbors, 
next-neighbors, next-to-next neighbors, and so on). Thus, a choice will have to be made to 
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see how many of them are useful in a given experiment or investigation. This certainly gives 
rise to an additional complication in the study of non-equilibrium system. 

Our approach allows us to associate affinity in a formal sense with all internal variables. 
This is how the classical non-equilibrium thermodynamics has been developed liQ- As 
observed by Landau and Lifshitz |2l], the use of internal variables in a modern way can 
be traced to Mandelstam and Leontovich jsG]; see also Pokrovski 37 1. Under the internal 
equilibrium assumption, Prigogine addresses the issue of internal variables (orientation of 
a molecule, deformation due to flow, elastic deformation, etc.) in Sect. 11, Chapter III 
of his classic book [12], or in Sect. 10.4 in the modern version (llj, and couples them to 
their "chemical potentials" or affinities. Indeed, Prigogine and Mazur were the first one 



to do this in their classic paper [38|; see also Coleman and Gurtin j39l]. The issue of the 
internal variables is also discussed in Sect. 6, Ch. 10 in Pokrovski [37| provides a very 
illuminating discussion of internal variables and their role in determining the internal energy. 
Thus, we will treat internal variables as additional thermodynamic extensive quantities or 
" observables" similar to the number of chemical species in chemical reactions that can be 
controlled by affinities or chemical potentials. More recently, the idea has also been visited 
by Bouchbinder and Langer j4|. 



VIII. THERMODYNAMICS OF A SIMPLE ROTATING BODY 



A. General Case 



We will find it convenient for later use to consider observing a body in different frames 
of reference; see also Appendices |X] and [Bl For concreteness, we consider the system E and 
assume that no internal variables and no other observables besides the energy, volume and 
number of particles are present; the latter can be added easily as we will discuss later. We 
will consider three special frames: the lab frame denoted by £, an intermediate frame X, 
with its axes parallel to those of and moving with respect to £ with a velocity V(t), and 
a frame C with its origin common with X and rotating with respect to it with an angular 
velocity f2(t). Let R(t) denote the location of the origins of X and C in the lab frame C at 
time t with R(t = 0) = 0. Let rc(t) denote the coordinate of a particle of S in the C frame, 
and vc(t) its velocity in this frame at time t. Its coordinate in the lab frame C is given 
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by 

= R + rc; (74) 

its velocity is given by Eq. (]A1|) . As shown in the Appendix [Aj the energy of the particle 
in the two frames are related as shown in Eq. (IA4p . Let us consider X to be the frame in 
which the center of mass of the body is at the origin, we will call it the center of mass frame 
for the body. Then, applying the above two relations to all the particles in the system and 
averaging over all allowed microstates [jj], which is carried out later in Sect. IVIII CI we 
obtain that the energy of the system in the three frames are related as shown in Eqs. (IBSp 
and flB8|) : 

Ec = Ec-^-M-n =Ex-M ■ n, (75a) 

p2 

Ex = Ec TT, (75b) 

where P and M are introduced in Eq. ( 1B2I) . We have not used the overbar to express the 
statistical averages as explained in Sect. IVIII CI but is implied. 
We first prove the following theorem: 

Theorem 5 The entropy of a system is the same in all three frames C,X, and C. 

Proof. To prove the theorem, we proceed as follows. Only for simplicity of the argument 
and presentation, we focus on a system with fixed V = V{t) and at some instant t. The 
extension to considering other extensive variables is trivial. Consider observing the system 
simultaneously at t in these frames. It is evident that corresponding to each pair r£, of 
the coordinates and momenta of a given particle in the lab frame C at this moment, there 
is a unique pair rj, pj and r^, p^ in the other two frames. This is true of all the particles. 
The collection of positions and momenta of all the particles defines a point in the phase 
space. In classical statistical mechanics, a microstate of the system is identified by a small 
volume element of size (27r^)^^ about a point in the phase space. Thus, corresponding to 
each microstate i (= ic-, ix, or ic) in one frame, there exists a unique microstate in the other 
two frames. The uniqueness of microstate-mapping ensures that their probabilities in the 
three frames are also equal: 

Pic = Pii = Pic- (76) 

Let us consider all the microstates of the same energy Ec in the C frame at time t, and 
let W{t) = W{Ec,t) denote their number and Pc{i) the set of their probabilities (not to 
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be confused with momenta pc, etc.). Because of the uniqueness of the mapping of these 
microstates noted above, not only the number of microstates in the three frames are the 
same at that instant 

W{t) = WiEct) = W{Ec,t) = W{Ei,t), 
but also the set of their probabilities 

Pcit) = Pci^) = Piit); 



however, their energies are different as given by Eqs. fl75l) fl75ap and (175b|) . This immediately 



shows that the entropies using the general Gibbs formulation in Eq. fl2Tl) are equal in the 
three frames: 

Sc{Ec,t) = Sx{Ej,t) = SciEct), (77) 

whether the system is in internal equilibrium or not. This proves the theorem. ■ 

It should be noted that the center-of-mass kinetic energy P^/2M is the same for all 
microstates in W{t). Similarly, it follows from Eq. (IB 10 1) that even M ■ f2 is the same for all 
microstates in W(t); see also the discussion leading to Eq. (]84b[) . Thus, the three energies 
only differ by some constants at each instant t. 



B. System under Internal Equilibrium 

We now specialize and assume the existence of the internal equilibrium, so that all mi- 
crostates are equally probable 

p,{t) = 1/W{t),z = 1,2,- ■■ ,W{t). (78) 

Hence, the three entropies are each equal to 

S{t) = In W{t). (79) 

It follows from Eq. fl77|) that there is no reason to use different subscripts to distinguish 
the entropies. Accordingly, we will use 5* to represent the entropies in different frames; their 
energy arguments will of course depend on the frame of reference. The arguments V and 
fl above are actually external parameters that are not extensive. We will show below that 
the entropies in the X and C frames are actually functions of extensive quantities P and M 
that are conjugate to V and fl, respectively; cf. Eq. (1271) . 
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C. Statistical Averaging over Allowed Microstates 



We now investigate the consequences of statistical averaging over microstates with non- 



zero probabihties 15|] and show that its consequences are the same as expressed in Eqs. fl75]) 
and (1771) . We first note that M in Eq. (lB2p depends on the coordinates and momenta of 
the particles, but this is not the case with P, even though both are extensive quantities. As 
Ec in Eq. flB5l) or fl75ap is for a microstate determined by the coordinates and momenta of 
the particles, we need to average it using microstate probabilities in Eq. (1781) . Averaging 
over various microstates relates the average energies in the two frame. We use an overbar, 
see Eqs. ( 12^ and to denote the average. We find that the same form also describes 
the desired relation between the average energies: 

Ec{t) = Ec{t) - ^ - M(t) ■ n{t) = Exit) -M{t).n{t), (80) 

where 

Exit) = Ecit) - ^. (81) 

The momentum P, of course, does not require any averaging as noted above. Eq. flHOj) is 
valid at each instance t. We can also take the statistical average of Eq. f lB3a|) to obtain 

Ec{t) = Ec{t)-^- ^m,r, ■ (v, x 17) - i^m,(17 x r,)^, (82) 

i j 

where the two sums are over all the particles in the system. Here rj and Vj are the instanta- 
neous position and velocity of the jth particle in a microstate with respect to the C frame; 
we have suppressed the subscript C from and Vj for the sake of notational simplicity. In 
the last equation, the third contribution is due to the relative motion of the particles with 
respect to the C frame. Indeed, the average of Eq. flB4p immediately yields 



M{t)-Q{t) = ■ (v, X 17) + J]m,(17 X r,)2. (83) 

j j 

The third contribution in Eq. fl82|) and the first contribution in Eq. fj83|) vanish when the 
system is in internal equilibrium because of the absence of any relative motion in that case; 
see Theorem [H 

Since Ej also does not depend on the velocity V, a similar averaging of Eqs. ( 1B9P and 
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(IBlOp gives us 

(if).,,..-' 

Comparing the above equations with the equations in the Appendix [B], we see that there is 
no reason to make a distinction between M(t), used in the proof above, and M(t) or the 
average energies and the energy used above in the proof. This justifies not using overbars 
to indicate statistical averages in Eq. (175|) . 

Since the entropy S(t) in Eq. fl7^ is fixed for fixed Ex, V, N, and il, we can express the 
above two derivatives at fixed 5* instead of fixed Ex- 

(if)....--'"' 

The above equation is similar to the well known result [16, Sect. 11] in equilibrium statistical 
mechanics that the statistical average of the derivatives of the the energy with respect to 
external parameters (V and Q) should be taken at constant entropy and other extensive 
quantities. We have extended this result to internal equilibrium now. Introducing the 
following standard derivatives 



nt), = -p{t) (86) 



defining the temperature and pressure of the system, we can write down the following dif- 
ferential identity 

dEc = T{t)dS{t) - P{t)dV{t) - M (t) ■ dn{t). (87) 

It should be noted that because of Eq. fl84al) . the average energy Ec(t) does not depend on 
the velocity of the frames X and C. Thus, there is no reason to keep V fixed in the various 
derivatives in Eqs. ( I8lll86|l . 
For 

Ex = Ec + M-n, 

we find that 

dEx = T{t)dS{t) - P{t)dV{t) + Vt{t) ■ rfM(t), (88) 
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which is an extension of the result given in Landau and Lifshitz [16|, Sect. 26] to the internal 
equilibrium. The point to note is that the entropy S{t) in the X frame is a function of the 
conjugate variable M(t) instead of f2(t). However, for 

Ec{t) = Ec{t) + ^ + M{t).n{t), 

we also find an additional contribution due to V: 

dEc{t) = T{t)dS{t) - P{t)dV{t) + V(t) • d-p{t)+n{t) ■ dM{t) (89) 

given in terms of all extensive quantities. The additional contribution due to the momentum 
differential dP{t) is due to the velocity of the system as a whole and is important to include 
in the lab frame. For example, such a contribution is needed to describe the flow of a 
superfiuid in which the normal and superfiuid components have different velocities so that 
the superfiuid cannot be considered at rest in any frame 2I, see Eq. (130.9)]. We will need to 
allow for this possibility when we extend our approach of nonequilibrium thermodynamics to 
inhomogeneous systems where different subsystems will undergo relative motion. It follows 
form Eq. (!89l) that the drift velocity of the center of mass of the system is given by 

if) --^w- (-) 

y'^) / S,V,N,M 

Similarly, the angular velocity is given by 

( dEcm 
\dM{t)J 



= n{t). (91) 

S,V,N,P 



We again observe that the entropy in the lab frame £ is a function of the extensive 
conjugate quantities P(t) and M(t) rather than the external parameters V and fi. 

From now on, we will not use the overbar to show statistical averages for the sake of 
notational simplicity. 

It is clear from Eq. ( IHTl) that we must treat Ec{t) as a function of S{t),V{t) and Q{t) 
for constant A^. Alternatively, we must treat S{t) as a function of Ec{t), V{t) and f2(t) : 

Scit)^SciEcit),Vit),nit),N), (92) 

which is identical to the functional dependence shown in Eq. fl77|) . except that we no longer 
have an explicit t-dependence because of internal equilibrium. The important point to 
observe is that the entropy is a function of not only the energy in the C frame, but is also a 
function of the angular velocity of the reference frame when rotation is involved. 
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D. Same Temperature and Pressure in Different Frames 



We now make an important observation. It follows from Eqs. ( 188|) and (!89|) that 

( W) ^,^,M ^ ( W) V^,^,P,M ^ ^^^^ 

obtained by differentiating with respect to S{t). Similar equations are obtained when we 
differentiate with respect to V{t). 



(dE^\ ^ fdEcit) 



= -m- (94) 

S,N,M \ V"; / S,N,V,M 



These equations are identical to the derivatives in Eq. (1861) and show that the internal 
temperature T(t) and the internal pressure P{t) are the same in the three frames. Moreover, 
it is the same entropy function S(t) that appears in Eq. (IHTj) also appears in Eqs. fl88l) and 
( 189|1 . In other words, the entropy is the same in all frames, except that the arguments are 
different. 



IX. THERMODYNAMICS POTENTIALS AND GIBBS FUNDAMENTAL RELA- 
TION FOR A HOMOGENEOUS SYSTEM WITH TRANSLATIONAL MOTION 

A. Thermodynamic Potentials for a System under Arbitrary Conditions 

1. Fixed number of particles N of the system S 

Before we discuss the inhomogeneous case, let us consider the homogeneous situation 
considered in I and revisited briefly in Sect. IVIl and extend it to the case when the system 
S moves as a whole with a linear momentum P. We still assume that Sq is at rest. Because 
of the linear momentum conservation, the linear momentum of the center of mass of S is 
—P. Thus, the centers of mass S and E are moving towards each other. For simplicity, we 
will assume the absence of overall intrinsic rotation for S and S individually. This can easily 
be incorporated as we do in the next section. Then we only need to consider the orbital 
angular momentum Lq of Sq in terms of the locations R and R of the centers of mass of S 
and S, respectively. It is clear that Lq always vanishes since the centers of mass of S and S 
are moving towards or away from each other so that P and R— R are colinear : 

R X P-R X P = 0. 
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The (average) internal energies of S and S in their center of mass frames (the C frame) 

are 

E' = E- PV2M, 
& = E- PV2M, 

while E and E denote their total energies in the lab frame C, respectively; see Eq. ( 1751) . 
However, because of the extreme large size of S, its mass M satisfies the inequality M » M, 
so that we can replace E^ by E without any appreciable error. The entropy of S is a 
function of the internal energy E^; however, this is not relevant for our argument here if we 
are only interested in identifying the appropriate thermodynamic potential for the system. 
The energy 

Eq = E + E 

of So remains constant in time. As discussed above, the additivity of energy is valid under 
the assumption that the interaction energy Eq^^^ (t) between the system and the medium is 
negligible. This ensures that the entropies are also additive. In the lab frame C, the entropies 
of the S and S are obtained by considering their entropies in respective rest frames C and 
C; they are S{E\ V, N, t) and S{E^, V, N, t), respectively; recall that we have set = for 
each of them. 

Using the fact that the medium is under internal equilibrium, we modify Eq. (1531) to 
reflect the dependence on internal energies to obtain 

So{Eo, Vo, iVo, t) = S{E\ V, iV, t) + S{E\ V, N) ^ S{E\ V, iV, t) + S{E, V, N). (95) 

We now expand and follow the steps in arriving at Eq. (I55l) : the steps are unaffected by the 
motion of S. We thus obtain 

S'o(t) -S- S{E\ V, N, t) - E{t)/To - PoV{t)/To 

in terms of the energy and volume of the system. We can now identify the Gibbs free energy 
and enthalpy in the lab frame C in terms of the energy of the system: 

G{t) = E{t) - ToS{t) + PoV{t), H{t) = E{t) + P,V{t)- (96) 

compare with the Gibbs free energy in Eq. (l52l) . Thus, the second law in terms of the Gibbs 
free energy remains unchanged and is given by Eq. (|58|) . 
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In the center of mass frame C of the system, the Gibbs free energy and the enthalpy of 
the system are given by 

G\t) = E\t) - ToSit) + PoV{t), H\t) = E\t) + P^V{t). (97) 

Note that the above functions depend on the internal energy and not the energy of the 
system S. But they are not useful in the lab frame in which the system is being observed. 
Thus, we conclude that the overall motion of the system does not change the enthalpy and 
the Gibbs free energy; we must use the appropriate energy in the frame of observation; the 
temperature and the pressure of the medium are not affected by the choice of the frame as 
noted near the end of Sect. IVIIII Similarly, the entropy of the system is unaffected by the 
choice of the frame as shown by Theorem |5l 

2. Fixed volume V of the system S 

Let us assume that instead of keeping fixed, we keep the volume of the system fixed. 
Then following the procedure given in Sect. W\\ we find find that the correct thermodynamic 
potential now in the two frames C and C are given by 

n{t) = E{t) - T^s{t) ~ ^Ji^N{t), 

n'{t) = E'{t) -ToS{t) ~ fioN{t), 

respectively. 

3. Fixed N and V of the system E 

Let us assume that we keep and V fixed. Then following the procedure given in Sect. 
IVH we find find that the correct thermodynamic potential now in the two frames C and C 
are given by the Helmholtz free energy 

F{t) = E{t)-ToS{t), 
F\t)=E\t)-ToS{t), 

respectively. 
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4- Extension to many state variables 



From now on, we will list energy, volume and particle number for any body separately 
and use X and Z to denote the rest of the observables and state variables, respectively. We 
will only fix the number of particles A^, but allow all other state variables to fluctuate. In 
this case, S{E\ V, N) in Eq. ( 195|) must be written as S{E\ V, N, Z) and its expansion in 
terms of small quantities gives 



S{E\ V, N, Z) ~ S{Eo, Vo, N, Zq) - | S 

oE 



Z(t). 



Here, |g corresponds evaluating the derivative at Eq, Vq, Nq, Zq (Xq and Iq) so that these 
derivatives are constant, independent of the properties of the system; see the discussion in 
deriving Eq. ( 15^ . Introducing the corresponding "chemical potential vector" /Xq for Xq 
and the "affinity vector" slq = Aq/Tq = (see Sect. IVIII) for Iq because of the internal 
equilibrium of E 



dS 



Mo dS 
To ' Ul 



Ao 

To 



we can identify a new thermodynamic potential G'^it) 
S{Eo,Vo,N,Zo)] : 



0, (98) 

-ro[S'o(^o,K,,iVo,Zo,t) - 



G'^(t) = E{t) - ToS{t) + PoV{t) + Mo ■ X(t) = G{t) + /Xq • X(t) 



(99) 



in the lab frame C As Eo,Vo, Nq, Xq and are constant, we have 



dlojt) 
dt 







because ao = 0. Thus, S{Eq, Vq, N, Zq) is a constant, and the second law tells us that 

dSo 1 dG^ 



dt 



To dt ' 



(100) 



as expected in any spontaneous process. In the C frame, we will instead have 



G'^(t) = E\t) - ToS{t) + PoV{t) + Mo ■ X(t) = G\t) + Mo • X(t). 



(101) 



The important point to note is that chemical potential vector /Xq the affinity vector Ao = 
are associated with the medium, just as Tq, Po are. The analogue of the thermodynamic 
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potential Q{t) and F{t) are 

n^(t) = n{t) + ■ x(t), n^{t) = n\t) + ■ x(t), 

= Fit) + ti, ■ X(t), F'^(t) = F'it) + Mo • X(t), 

respectively. Again, it follows from the second law that 

dQ^ dF^ , , 

as expected in any spontaneous process. 

B. Gibbs Fundamental Relation for a System under Internal Equilibrium 

1. No extra state variables Z(t) 

We will first assume that there are no internal variables, but the system S satisfies the 
condition of internal equilibrium so that S{E^, V, N) no longer depends explicitly on time. 
Then, we can identify the temperature, pressure, and the chemical potential of the system 
by 

dS\ 1 



dE' J T{t) ' 



(103a) 



(103c) 



dN J T{t) 

This allows us to write down the Gibbs fundamental relation for constant as 

T{t)dS = dE\t) + P{t)dV{t) - +i^{t)dN{t), 

which can be rearranged to write down the first law of thermodynamics 

dE'{t) = T{t)dS - P{t)dV{t) + ij{t)dN{t). (104) 

We now turn to Eq. ( IHOl) . In terms of the momentum of the center of mass frame C of the 
system, we have 

dE' = dE-Y ■ dP; (105) 
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recall that according to our assumption, the system has no intrinsic angular momentum. 
This allows us to use Eq. fll04p to write down the differential form 



dE = T{t)dS + V ■ rfP - P{t)dV{t) + n{t)dN{t) 



(106) 



for the first law of thermodynamics in terms of the energy E{t) rather than the internal 
energy E\t). This has some important consequences and will be extremely useful in the 
following. The first consequence is that it allows us to think of S{E\ V, N) as a function of 
four variables S{E, P, V, N): 



This equation is the Gibbs fundamental relation relating the entropy with E(t) rather than 
E\t). The second consequence, which follows from the Gibbs fundamental relation is that 



The drift velocity V (of the center of mass) of the system is given a thermodynamic inter- 
pretation in terms of the derivative in Eq. ( llOTbj) at fixed E, V and A^. For the case when 
the number of particles is held fixed, dN{t) = and the last term in dE\t), dE{t) and dS{t) 
will be absent. For fixed V , the third term in dE\t), dE{t) and dS{t) will be absent, and so 
on. 

2. Inclusion of state variables Z{t) 

In the presence of internal variables X and I, the extension of the above relations is quite 
obvious. Introducing the instantaneous chemical potential vector ^t(t) associated with X 
and the affinity vector A associated with I using 



T{t)dS = dE~Y ■dP + P{t)dV{t) - ii{t)dN{t). 




(107d) 



(107b) 



(107a) 



(107c) 




(108) 
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we generalize the entropy differential to 



T{t)dS = dE\t) + P{t)dV{t) - ii{t)dN{t) - ^i{t) ■ dX{t) + A(t) • dl{t) 



(109a) 



dE{t) - V(t) ■ dP{t) + P{t)dV{t) - ii{t)dN{t) - nit) ■ X(t) + A(t) ■ dl{t). 



(109b) 



X. INHOMOGENEOUS SYSTEM WITH RELATIVE MOTION 
A. Subsystems undergoing Relative Motion 

We now consider the isolated system to be stationary so that it has no linear and angular 
momenta. However, S and S may have relative motion so that they each may possess linear 
and angular momenta that individually must cancel out: 



Moreover, we will now treat E as inhomogeneous and assume that it can be decomposed 
into a collection of a large number No, of subsystems a^i k = 1,2, ■ ■ ■ , N^, which may be in 
different macrostates to allow for inhomogeneity and for relative motion between different 
subsystems and within each subsystem. Each subsystem is still macro scopically large so 
that we can not only introduce a legitimate entropy function Sk for the macrostate M.k via 
Gibbs's formulation 



where au denotes one of the allowed microstates of the subsystem corresponding to 
the macrostate M.^ characterized by observables e\,nk,Vk and z^, but we also have these 
entropies satisfy the additive property 



which requires their quasi-independence as discussed in Sect. |V] at each instant t. Using 
the entropy s^, we can introduce the appropriate thermodynamic functions, but care must 
be exercised in identifying these functions in the lab frame C, where the experiments are 
done. This is because the energies depend on the frame of reference, which will result in 
different values of the energies and thermodynamic potentials in different frames, such as 



P + P = 0, M + M = 0. 






(Ill) 



k=l 
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the lab frame L and the rotating frame frame attached to the center of mass of 0"^, which 
is translating with a linear velocity Vfc(t) and rotating with an angular velocity i^fc(^)- 

To make further progress, we will make another assumption later as we did in I that each 
subsystem is in internal equilibrium. This occurs when all microstates contributing to the 
entropy are equiprobable 

Pa,{t) = l/Wk, a = 1,2,,- ■■Wk{t); 

here Wk{t) represents the number of microstates of the subsystem at time t. Under the 
equiprobability assumption or the internal equilibrium assumption, 

Skit) = In Wk{t), 

which is what one would obtain by applying the Boltzmann formulation of the entropy 
[isl . It also follows from Theorem [1] that the entire subsystem is uniformly translating with 
a linear velocity Vk{t) and rotating with an angular velocity ^lk{t) so that the system in 
internal equilibrium is stationary in the frame Ck- 

B. System under Arbitrary Conditions 

The internal energy for each subsystem ak is related to its energy Ck in the lab frame C 

el = Ck - pl/2mk - nik ■ flk, (112) 
see Eq. ( !80l) . Alternatively, we can use Eq. ( !82l) to express this relation as 

4 = ek- Pfc/2mfc - ' (^i^ ^ ~ \ ^ ^^^2) 

jk jk 

where each sum is over Uk particles in the subsystem k. The third term in the last equation 
vanishes when the subsystem is in internal equilibrium as commented earlier. We also have 
the additivity laws 

l^ = ^t;fc,Ar = ^nfc,X = ^Xfc,I = ^ifc, (114a) 
£^ = 5^efc = ^(e^ + p^/2mfc + mfc-fifc), (114b) 
P = Y^Pk,M = J2i^k + rkXpk), (114c) 



59 



at each instant t. The angular momentum 1^ = x is the orbital angular momentum of ak 
with Tfc, pfc representing the location and momentum of the center of mass of ak, respectively, 
and should not be confused with its intrinsic angular momentum rrifc introduced in Eq. (lB2p . 
We should emphasize that the additivity of the energy requires that the interaction energy 
between subsystems be negligible. This condition is necessary for the entropy to be additive 
as discussed in Sect. |Vl As a consequence, we do not have the contribution analogous to 



ilj{r)dVmEq. (g]) 
thermodynamics 



. This distinguishes our approach with that taken in local non-equilibrium 



J3- 



We keep fixed for simplicity so that is also fixed, and allow e\, and to vary in 
time. We assume as above that S has well-defined field variables (Tq, Pq, etc.) which do not 
change with time. This is ensured by assuming internal equilibrium for the medium. This 
is the only assumption we make regarding the isolated system Eq for which we have 

Eo = E + E,Vo = V + V,No = N + N,Xo = X + X,lo = I + I. 

Because the medium is taken to be in internal equilibrium, its energy is related to its internal 
energy according to Eq. ( IHOl) 



E' = E- P^/2M -Mn = E- P^/2M - ^1^/21, (115) 

where according to Eq. flllOp . we have taken M = — M, and where / is the moment of 
inertia of the medium about its axis of rotation. The axis of rotation must be one of its 
principal axis of rotation; see the comment after the proof of Theorem [1] in Sect. ??. The 
contribution coming from the internal motion, which is similar to the third contribution in 
Eq. (11 13p applied to the medium vanishes because of its uniform translation and rotation 
following Theorem [1] applied to the medium. The contribution similar to the last term in 
Eq. (I113p is the standard rotational kinetic energy of the medium treated as a rigid body. 
The angular momentum M is given by 

Assuming that the motions are finite, we conclude that must be finite. Therefore, for 
an extremely large medium, O and V must be extremely small, which ensure that the last 
two terms in Eq. (IllSp are extremely small. This allows us to approximate 

E c^E' (116) 
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without any appreciable error. As far as the system is concerned, the relationship between 
its internal energy and the energy is still given by Eq. (IB5|) 

(117) 



E' = E- pV2m - m ■ rj, 



except that the system may not be uniformly translating and rotating about any of its 
principal axis of inertia. We take and as constants, but allow for E, V, and Z (X and 
I) to change with time. Now, we use the entropy additivity 



and expand 



^o(^i, A^o, Vo, Zo, t) = SiE\ n,N, V, Z, t) + SiE\ n,N, V, Z) 



dS 















SV 





dZ 



S{E\ n,N, V, Z) = SiEo, n,N, Vo, Zo)-E-^ 

oE 

= SiEo, h,N, Vo, Zo) - ^ - - X 

-fo -fo -fo 



Mo 



by treating E, V, X and I as small parameters; here, |q corresponds to evaluating the deriva- 
tives at Eo, fi,Vo, Ao, Zo (Xq and lo). We have also used the definitions of the field variables 
of the medium in terms of the derivatives of the entropy S in Eqs. (jSlj) and ( l98l) . and set 
Ao = as established in Sect. I VI II We thus finally obtain 



^o(^o, A^o, Vo, Zo, t) - S{Eon,N, Vo, Zo) = S-^-V^ 

~~%~2MTo' 



G Mo 
~-~To~^'To 
M ■ n X /Xo 



where 



G = E-ToS + PoV, G' = E'- TqS + PoV. 



Observe that the expansion of S{E\ fl,N, V, Z) does not require any knowledge of what 
is happening with the system S. For example, we have not assumed internal equilibrium 
in the system or any of its subsystems. Since S{Eo, N, Vo, Zo) is a constant as established 
earlier, the second law tells us that 

1 dG^ 



dSo 



> 0, 



(118) 



dt To dt 

the desired result for the system. The inhomogeneity of the system has no relevance in the 
above conclusion. 
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Let us now assume that X = and introduce 



gk = ek- %Sk + PoVk, gl = 4 - T^Sk + P^Vk 



(119) 



for the subsystem (Jk- It is clear that 



G = Y,9k.G' = Y,9k. 



k k 

Because of the quasi-independence of various subsystems, Eq. flllSp immediately leads to 

dgk 



dt 



< 0. 



(120) 



Thus, gk(t) can be identified as the Gibbs free energy of the subsystem ak in the lab frame 
C. Comparing this definition with the definition in Eq, ([6]) used in the conventional non- 
equilibrium thermodynamics, we see that the discrepancy in the two Gibbs free energy has 
not disappeared by taking into account the inhomogeneity inherent in the system. 

It may be argued that the above identification of gk in Eq. f lllQp is based on considering 
the entire system in contact with the medium. One can alternatively consider a particular 
subsystem a of the system in contact with the medium and the remaining subsystems. 
However, a little bit of reflection shows that this will not affect the behavior of gk for the 
simple reason that the remaining subsystems still form a very small part of the isolated 
system so that they cannot affect the internal equilibrium of the medium. To see this more 
clearly, let us introduce a new medium S' consisting of S and the remaining subsystems. 
Only for the sake of simplicity, we do not consider any additional state variables. The 
argument can be easily extended to include them. Let S' denote the entropy and E'-' ~ E', 
V and A^' the internal energy, volume and the number of particles of S', the latter of which 
is kept fixed. Expanding this entropy in terms of the small quantities e' and v of the chosen 
subsystem requires calculating the derivatives 



dS' 



dE' 



and 



dS' 



Eo,Vo 



dV 



Eo,Vo 



Because of the small size of the system relative to the original medium E, these derivatives 
are not different from 



dS 
dE 



, dS 
and 

dV 
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respectively. Thus, using e = Eq — E for the energy of the subsystem a^, we find that 



S'{E", n,N', V) = S{Eo, n,N, Vo) - e 



dS' 



dE' 



dS' 



Eo,Vo 



dV 



Eo,Vo 



SiEo,N,Vo)-^-v^. 



Using this in ^o(^0' ^o, Vo, t) = s{e\n,n,v,t) + S'{E'' ,n,N' ,V'), where O is the angular 
velocity of the C frame of the selected subsystem, allows us to identify 



g = e-TQS + PqV 



(121) 



as the Gibbs free energy of the particular subsystem a. It now follows from Eq. (|2]) that 

a property we expect from the Gibbs free energy of a system. Incidentally, this also provides 
an independent justification of the inequality in Eq. f ll20p . 
In terms of g],, we immediately have 



which is expected in view of the sum rule in Eqs. f lll4bp and f lll7p . 



C. System under Internal Equilibrium 

The above derivation only uses the second law, and the assumption that the medium sat- 
isfies the condition of internal equilibrium. The situation with the first law is very different. 
In general, the differential ds^ of the entropy of the subsystem ak is given by 

dsk i dsk dsk dsk ,„ dsk dsk 

dsk = T^—de^ + T^dvk + 7^ — duk + ■ dllk + dzk{t) + -Trrdt] 

ock ovk orik oilk 07.k ot 

cf. Eqs. (177|) and (!92l) . If we now assume internal equilibrium so that Sk does not have an 

explicit t-dependence, we can get rid of the last term above. In this case only, the remaining 

derivatives identify the field variables 1/Tk, Pk/Tk, —fJ^k/Tk etc. of the subsystem ak : 

dsk ^ 1 dsk ^ Pfc(t) dsk ^ Atfc(t) dsk ^ mfc(t) dsk ^ ^i^^\t) dsk ^ A(^')(t) 

dck ~ Tk{t) ' dvk ~ Tk{t) ' duk ~ Tk{t) d^k ~ Tk{t) ' d^k ~ Tk{t) ' dik ~ Tk{t) ' 

(122) 

These derivatives then allows us to obtain the Gibbs fundamental relation 
de[{t) = Tk{t)dsk-Pk{t)dvk + fik{t)dnk-uik{t)-dnk + tJ'^''\t)-d^k{t)-A'^''\t)-dik{^^ (123) 
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Using 

dek{t) = d^it) + Vfc(t) ■ dpk{t) + d[nkit) ■ mfc(t)], 
we find the first law of thermodynamics can be expressed in terms of the energy as 

dek = Tk{t)dsk{t) + Vfc(t) ■ dpfe(t) + ^^fc(t) ■ rfmfc(t) - Pk{t)dvk{t) (124) 
+ lik{t)dnk{t) + ^W(t) ■ dxfc(t) - A('^)(t) ■ difc(t); 

compare with Eq. (!89|) . This allows us to think of the entropy Sfc((e)., Jlfc, n^, f^, z^, t) as a 
function 

Sfc(efc,Pfc,nifc,nfc,Wfc,Zfc,)f:) 

so that the drift and the angular velocities in internal equilibrium are given by 

Vfe(t) _ dskjt) Clkit) _ dskjt) 

T,{t) 9pfe(t)' T,{t) dm,{ty ^ ^ 

However, different subsystems will undergo relative motions with respect to each other as 
Vfc(t) and Ofc(t) are different for them. In addition, their temperatures Tk{t) and pressures 
Pk{t) are also different for each other, and so are fi^'^^t) and A^'^^t). Thus, there would be 
viscous dissipation and, consequently, entropy production as they come to equilibrium with 
each other and with the medium. We now turn to this issue. 



XI. REVERSIBLE AND IRREVERSIBLE CONTRIBUTIONS 
A. General Considerations: Bodies in a Medium 

It is a well-known fact [IG,, Sect. 10] that in equilibrium, the system E has a uniform 
translational motion as a whole with a constant velocity Vq and a uniform rotation of 
the whole system with a constant angular velocity fio- Thus, there cannot be relative 
motions between different subsystems in equilibrium. Thus, no internal macroscopic motion 
is possible in equilibrium. In equilibrium, the coefficients of the differential quantities in Eq. 
(11241) take their constant values of the medium. Writing 

dSkit) = deSkit) + diSk{t) (126) 

as a sum of the change in the entropy deSk{t) due to reversible exchange with the medium 
and the production of the entropy diSk{t) > due to irreversible processes within the system, 
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we have 



deSk = ^[dek{t) - Vo(t) ■ dpkit) - flo ■ dnikit) + Podvk{t) - fiodukit) - /Xg ■ '^Xfe(t)] 



and 

diSk = F 
+ F 
>0, 
where 



[Tkit) 
Tkit) 



dskit) + F 
dukit) + F 



Vfc(t) 



Tkit) 



Tkit) 



dpkit) + F 
■d^kit) + F 



^kit) 
Tkit) 

Tkit) 



■ dnikit) + F 

■ dikit) 



Pkit) 
Tkit) 



dvkit) 

(127) 
(128) 

(129) 



F^=F [w] = wit) - Wcq 

is the difference of the conjugate field wit) at time t and its value in equilibrium, i.e. as 
t —7- oo, and represents the thermodynamic force associated with the conjugate extensive 
quantity ^(t); see also Eq. (ITS!) . For the internal variables, the equilibrium value of Aq is 
zero. According to the second law of thermodynamics, not only d^Sk > 0, but each term in 
Eq. f ll27p is non-negative: 

F [w] dz > 0, 

where z and w form a conjugate pair. In terms of these pairs, we can express the two parts 
of the entropy as follows: 



dsk 



Wfc ■ dzk, deSk = wofe ■ dzk, diSk = F [w^] • dzk, 



(130) 



which is the general form of the entropy differenmtial and its two components. 

Let us now turn back to the current case under investigation. We can express the gener- 
alized Gibbs fundamental relation as 



dck = deCkit) + dfikit), 



where 



and 



d^Ck = Todskit) + Vo(t) ■ dpkit) + flo ■ dnikit) - Podvkit) + fiodukit) + /Xg ■ cixfc(t), 

diCk = F [Tkit)] dskit) + F [^rkit)] ■ dpkit) + F [^kit)] ■ dnikit) - F [Pkit)] dvkit) 
+ F [fikit)] dukit) + F [/xW] ■ d^kit) - F [A^^^] ■ dikit). 
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The generalized form for the Gibbs fundamental relation is evidently 

dtk = Tkdsk - Wfc ■ (iz', deCk = T^dsk - Wofe ■ rfz', d^Ck = F [Tk] dsk - F [Wfc] ■ dz', 

where z' represents all state variables except the energy e^, and where = Tw and 
Wo= ToWq. 

It is easy to see that performing the Legendre transform to obtain the Gibbs free energy 
in Eq. f lllQp only affects the form of deSk, but leaves diSk unaffected. Thus, it is easy to see 
that 

dgk(t) = degkit) + digk{t), 

where 

dc9k = -Sk{t)dTo + Vo(t) ■ rfpfc(t) + ■ dmk{t) + Vk{t)dPQ + fXodnk{t) + ■ rfxfc(t), 

and 

digk = diCk- 

The general form of dggk is 

degk = -Sk{t)dTo + Vk{t)dPo - Wfc ■ dz" , 
where z" represents all state variables except the energy Cfc and the volume Vk- 

B. Bodies forming a Finite Isolated System without a Medium 

So far, we have considered a system or a collection of subsystems in a very large medium 
S so that its field variables are held fixed at Yq. We now consider a collection of bodies of 
comparable sizes forming an isolated system Sq. In this case, we cannot treat any of the 
bodies as a (macroscopically extensively large) medium with a fixed field at Yq. As the 
collection strives towards equilibrium, their field variables continue to change in time. This 
causes a problem in identifying reversible contributions to various quantities. To solve this 
problem, we discuss below a simple case; the generalization to more complex situation will 
be obvious. 
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1. Two Bodies 

Let us consider the simplest possible case of two comparable bodies 1 and 2 in internal 
equilibrium at temperature Ti and T2 > Ti, respectively, that are brought in thermal contact 
at time t = 0. We will simplify the discussion by assuming that all other extensive observables 
besides the energy are held fixed. The case of two bodies in the shape of rectangular cuboid 
along the x-axis are shown schematically in Fig. |3]^a), with the rectangular interface of area 
A lying in the ?/2;-plane. Let dQ = dE be the infinitesimal heat given to the body 1 by the 
body 2. As the amount of heat is infinitesimal, it does not alter the temperatures of the 
bodies in any significant way so that the entropy of the isolated system Eq changes by 



at t = 0.; compare with Eq. ( I29l) . This expression will not be correct as time goes on as 
the temperature continues to change. The general expression for the irreversible entropy 
generation is 



where dQ{t) is the infinitesimal heat given to the body 1 by the body 2 between t and 
t + dt, and Ti{t), T2{t) are the instantaneous internal temperatures of the two bodies. It is 
clear from the discussion following Eq. (1291) that the irreversible entropy generation in Eq. 
(11311) is the sum of entropy changes in the two bodies. Thus, Eq. (11311) is exact for the 
isolated system consisting of the two bodies. However, the point to remember is that it is 
the irreversible entropy generation in the two bodies, and says nothing about the irreversible 
entropy generation within each body. 

If the system above were not isolated, then Eq. (I13ip would still give the entropy change 
due to the direct flow of heat between the two bodies, but it not represent the irreversible 
entropy generation due to this heat flow. Moreover, there would also be changes in the 
entropy of the system and each of the bodies due to heat exchange with the medium or other 
bodies. These entropy changes will also have their own irreversible entropy generations. The 
presence of the medium at the equilibrium temperature Tq of the isolated system is discussed 
below. The situation when To is not the equilibrium temperature of the isolated system is 
very different, as discussed later; see the discussion after Eq. ( 11321) . 

What can we say about the irreversible entropy generation within each body? This 





(131) 
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1: T,{t) 



T 
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2: 7^2(0 
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(b) 
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FIG. 3: We show in Fig. (a) a simplified situation in which two bodies 1 and 2 in thermal contact 
are alligned along the rc-axis. Initially, the two bodies are in equilibrium at temperatures Ti and 
T2 > Ti. Their contact interface is taken to be a plane orthogonal to the x-axis as shown in the 
figure, but the discussion is valid for any shape of the interface. In reality, however, any interface is 
an interface region of some very small width Ax over which the temperature continuously changes 
between Ti(t) and T2{t) so that there is a very narrow region of width dx « Ax around a point E 
on the rc-axis, where the temperature is exactly Tq, the equilibrium temperature of the two bodies, 
as shown in Fig. (b) for the simple case shown in Fig. (a). 



question is, to the best of our knowledge, is not answered within the local thermodynamics 
approach. To answer this question in our approach, we proceed as follows. We know that 
the entropy generation in each body must vanish when the bodies come to equilibrium. To 
obtain the desired result, we introduce the common temperature Tq, when the two bodies 
come to equilibrium. We now discuss two alternative approaches to determine the individual 
entropy generation on the way to prove Theorem [2l 
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a. Introduction of a Medium at constant Tq Let us imagine that we insert the two 
bodies in an extensively large medium S kept at a fixed temperature Tq, with all three 
bodies now forming a new isolated system Sq. The situation is schematically shown in Fig. 
m There is no direct contact between the two bodies, as shown. Despite this, the heat loss 
dQ(t) > by the body 2 to S is completely transferred to body 1. An alternative is to 
insert the medium between the two bodies along the x-axis but not surrounding them from 
all sides with the same effect of transferring the entire heat dQ{t) from body 2 to body 1. 
The width of the medium along the x-axis may be infinitesimally small of order dx, but 
must have a macroscopically large cross-sectional area in the yz-plane to ensure its constant 
temperature Tq at all times. In either case, the entropy of the medium does not change so 
that dS = 0, deS = 0, and d^S = 0. Therefore, the irreversible entropy generation in the new 
isolated system Sq is identical to the irreversible entropy generation in Sq. This artificial 
introduction of S can now be exploited to obtain the irreversible entropy generation in each 
smaller body. The method of calculation above can be easily applied to this case, see Eq. 
([I27D, to yield 

dMt) = dQit) - 1^ , dMt) = dQ{t) (^^ - ^) , (132) 

where dQ{t) > is the heat addede to body 1 or given out by body 2. By comparing with 
Fig. |2l we note that dQ(t) = dQi{t) + dQ'{t) = dQiit) + dQ'{t), giving 

dQiit) = dQ2{t). 

It is evident that these entropy productions vanish in equilibrium, as required by the notion 
of equilibrium between the two bodies. If we had introduced a medium held at a constant 
temperature Tq 7^ Tq, then the equilibrium will occur at Tq, and not at Tq, and we would 
end up with a different final state of Sq than that of Sq. Thus, the medium must have the 
constant temperature Tq. 

The situation will be very different if the medium is taken to be at a themperature 
Tq 7^ Tq. In this case, the heat given to body 1 is different from the heat given out by body 
2. Comparing with Fig. [2], we note for the total amount of heat dQ2^tot{t) given out by 
2 and the total amount of heat dQi tot(t) given to 1 that dQi^totit) = dQi(t) + dQ'(t) 7^ 
dQ2,tot{t) = dQ2{t) + dQ'{t), giving 

dQi{t) 7^ dQ2it). 
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FIG. 4: A schematic representation of a new isolated system Sq consisting of the two bodies 1 and 
2, not in thermal conatct, surrounded by an extensively large medium S, which is so large that 
the presence of 1 and 2 does not alter its internal temperature, which therefore remains constant. 
We choose the medium, as exlained in the text, to be at the temperature Tq. 



The irreversible entropy generation in the two bodies are now givem by 
so that diS{t) = diSiit) + diSi{t) is given by 

^■S(0 ^ ^Q'm (5^ - + - ^) ^QM - 5^) ■ (133) 

Notice that for > T2{t) > Ti{t), dQi{t) > 0, and dQ2{t) < 0. Similarly, for < Ti(t), 
dQi{t) < 0, and dQ2{t) < 0. Thus, each of the last two irreversible entropy generation 
contributions above is non-negative, as expected. We will see below that the last two con- 
tributions are absent in the local theory. 

The above approach can now be extended to many bodies at different initial temperatures 
Tfc. We assume that the medium S surrounds all the bodies so that there is no direct contact 
between the bodies. If Tq still denotes the final temperature of all the bodies, treated as 
an isolated system Eq, then the medium S must be chosen to be also at the same constant 
temperature Tq. It is obvious that 

d-Mt) = dQ kit) (^^^ - , (134) 
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where dQk{t) denotes the infinitesimal heat added to the fc-th body between t and t + dt, and 
Tfc(t) is its instantaneous internal temperature; compare with Eq. fll27p . This is a general- 
ization of Eq. (IT7|) to many bodies, and exploits the absence of causality inherent in deSk(t), 
as exemplified by Eq. f ITB]) . The irreversible entropy generation in the isolated system Eq 
consisting only of the bodies (without the medium) is the sum of all these contributions: 

k k J-kyl^) k 

as it must be for an isolated system; Obviusly, diSo{t) does not depend on the final temper- 
ature Tq. Here, we have used the fact that 

ZdQk{t) = 

k 

because of the isolation of Eg and Eq. It is quite clear that the discussion is easily extended 
to include other extensive variables which results in the expression for diSk{t), which is 
identical in form to the expression diSk{t)in Eq. f ll27p . 

One may feel uneasy that the introduction of the fictitious medium E has changed the 
problem. To see that this is not the case, we reemphasize that its introduction does not 
affect the heat dQk{t). Since, it is the heat dQk{t) that determines the entropy generation, 
and since the equilibrium state of each body is the correct equilibrium state, the result must 
be correct. To offer an even stronger argument, we now provide an alternative proof without 
the introduction of the medium so that we can feel comfortable in associating a medium in 
the case when finite-size bodies form an isolated system. The importance of this trick is 
that it allows us to use all the results we have found by the use of a medium. 

b. Without any Medium Let us revert to the simple case of two bodies in thermal 
contact. In reality, the interface or the contact region between the two bodies is not a sharp 
boundary with a discontinuity in temperature; rather, it is a narrow region of width Aa; over 
which the temperature rapidly changes from Ti(t) to T2(t), as shown in Fig. |3]^b). Thus, 
there exist a point E over this region where the temperature is precisely Tq. It is possible that 
the location of the point E along the x-axis changes with time, but this point is irrelevant. 
The relevant point is that the temperature in a very narrow width dx « Ax around and 
including this point will remain constant in time. If we take the narrow neighborhood of a 
point whose temperature is different from Tq, its temperature will eventually change to Tq, 
as equilibrium is achieved. Thus, the temperature of the narrow region around this point 
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will not be constant in time. We take the region where the temperature is greater than Tq 
along with the half- width dx/2 around E on the side of 2 to be the part of the body 2, 
while the region with the temperature less than Tq along with the half-width dx/2 around 
E on the side of 1 to be the part of the body 1. From the quasi- independence of the two 
bodies, it is clear that the inclusion of these regions will not significantly affect the internal 
temperatures Ti{t) to T2(t). The heat lost dQ{t) by the body 2 is transferred to the body 1 
at a constant temperature Tq. Thus, while 

^ ' Ti(t) ' 

where Ti{t) is the internal temperature of the body 1, 

Jo 

SO that 

4s,w = <;w)(5^-^). 

as discovered above; see Eq. fll32p . In a similar fashion, we obtain 

dMt) = dQ{t) ^ 



To m) 

in accordance with Eq. fll32p . 

The above discussion justifies the use of Tq as a temperature in the interfacial region. 
However, the important point is that the determination of deS{t) for any body requires the 
use of the equilibrium temperature of the body in accordance with Eq. ( TT6l) . Thus, the 
discussion is equally valid for any number of subsystems in the system. 

c. Comaprison with Local Thermodynamics The same result as in E g. (11311) is also 



obtained in the local thermodynamics, as can be easily seen; see for example [12|, Eq. (3.14)]. 
In the limit in which the interfacial region between the two subsystems is infinitesimal in 
width along the x-direction, we can obtain the continuum analog of the irreversible entropy 
generation between the two neighboring regions. Denoting Ti(t) by T{x,t) and T2{t) by 
T(x + dx, t) ~ T{x, t) + (dT/dx) dx, we have 

dSo{t) = d^Soit) - dQit) [d{l/T)/dx\ dx 

for conduction. Dividing and multiplying by the cross-sectional area A of the interface, 
and using Adx as the volume of the interfacial region, we have for the rate a of entropy 
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production per unit interfacial volume 

a{x,t)=q{t)[d{l/T)/dx\, (135) 

where q{t) = Q(t)/A denotes the heat flux. This expression (in three dimensions, the result 
will contain d{l/T)) in this limit is a standard result for the entropy production due to heat 
conduction in local thermodynamics. 

It should be stressed that our derivation of a{x,t) above is independent of the how long 
the two subsystems are along the x-axis, but assumes implicitly that T{x,t) has a Taylor 
series expansion. It is also obvious that a{x,t) must be zero outside the interfacial region. 
Thus, a{x,t) should be correctly identifled as the rate of entropy per unit volume in the 
interfacial region. Therefore, as it follows from the discussion immediately following Eq. 
f ll3ip . the derivation says nothing about how much of the irreversible entropy is generated 
within each body. The issue is avoided in local thermodynamics by assuming that the entire 
volume of the system is composed only of such interfacial regions. This is contrary to the 
basic postulate of local equilibrium according to which each local region has a well-deflned 
temperature T{t), while the interfacial regions have non-zero gradients of the temperature. 

We also observe that the form of a{x,t) in Eq. fll35p is only valid for the case when 
Eq. f ll3ip is valid. It is not valid for the case covered in Eq. fll33p . This is the case when 
the equilibrium temperature Tq of the system is different from the equilibrium temperature 
To of the two bodies under consideration if treated as isolated bodies. Thus, the above 
expression a{x, t) will not be valid if our system consists of a large number of bodies so that 
the irreversible entropy generation between any two bodies in contact will be given by Eq. 
f ll33p . In this case, the expresasion explicitly contains the equilibrium temperature of the 
system in the last two terms, which is not the case with a{x,t) in Eq. fll35p . thus verifying 
the earlier statment made after Eq. fll33p . 

2. Several Bodies 

Let us now consider a simple extension of the case shown in Fig. [3]^a) in which there are 
several bodies 1, 2, 3, ■ ■ ■ ,B in thermal contact along the a;-axis forming an isolated body Sq, 
with their temperatures in an increasing order: Ti < T2 < ■ ■ ■ < Tk < T^+i < ■ ■ ■ < Tb- Let 
Tk-i < Tq < Tk again denote the equilibrium temperature Tq of the isolated body Sq, so 
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that there exist a point E over their interface region where the temperature is precisely Tq. 
Let us consider the infinitesimal heat dQi{t) received by 1 to be the part of the heat dQi{t) 
that was given out by the body k and was transferred unaltered via k — 1, k — 2, ■ ■ - 3,2 to 
1. Any heat transfer through the interface region at constant To is isothermal. Therefore, 
the reversible entropy change due to dQi(t) is precisely d^Si = —dQiit) /Tq. The entropy 
change of 1 is dSi = dQi{t)/Ti{t) so that it immediately follows that 



for the body B, where dQ-Q{t) is the heat rejected by B, which was transferred isothermally 
to the body k unaltetred. For the body 2, which receives an infinitesimal heat dQ2{t), we 
similarly find 



and so on. All these results can be easily obtained by inserting a medium held at fixed 
temperature Tq at the interface between each connsecutive pair of bodies. Thus, we conclude 
that we can consider an isolated inhomogeneous body as a body embedded in a medium 
without affecting any of the consequences. This then proves Theorem [2j 

C. General Discussion 

We are now in a position to provide a proof of Theorem [2J We consider the reversible 
entropy change d^s^ in a given subsystem due to the p-th conjugate field W^{t) {Y{t) or 
A(t)) due to the change dz^. It is given by 




as above. Similarly, we find 






while 



ds" 



W{t) 



dzP, 



Tit) 




in terms of the quantity F [w] in Eq. (I129p . This thus proves the theorem. 
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XII. DISCUSSION AND CONCLUSIONS 



We have proposed a scheme to extend our previous work in I to an inhomogeneous 
system in which subsystems may undergo relative translational and rotational motions with 
respect to each other to contribute to irreversibility. Another source of irreversibility is 
the temperature variation, which was considered in I. Each subsystem or the medium is 
identified by a set of extensive state variables Z(t), and some constant parameters C, some 
of which may be external parameters and need not be extensive. Examples of C may be the 
number of particles in the system that characterize the system, or the translational and the 
angular velocities of the frame of reference in which the observations are made. Examples of 
the state variables are mechanical quantities such as energy, volume, etc. and the internal 
variables are the translational and angular momenta of the various subsystems, etc. 



A. Quasi-independence and Additivity of Entropy and Energy 

As our approach starts from the second law, the law of increase of entropy, the entropy 
as a state function plays the most important role in our approach. Accordingly, we need to 
ensure that this quantity can always be determined by or at least formally defined in terms 
of some basic quantities pertaining to the macrostate of the system of interest. We use the 
Gibbs formulation of the entropy such as that in Eq. ( I2T]) . which is applicable in all cases as 



discussed in a recent review 



15j. This formulation is well-established for an isolated system 



151], but we show in Sect. |V]that it is also applicable to an open system even when it is not 



in equilibrium with its surrounding medium. This formulation of entropy for open system is 



a well known result in equilibrium thermodynamics [16| and our demonstration generalizes 
this result to non-equilibrium systems so that this entropy can be used as the central quantity 
in developing a non-equilibrium thermodynamics with the second law as the starting point. 
It is required that the open system be quasi-independent of the medium. We have shown in 
Sect. |V]that this quasi-independence is a very important property, which is required for the 
additivity of the entropy and of energy. For example, quasi-independence of the system and 
the medium ensures that their energy of interaction is negligible so that dropping it makes 
their energies additive, as seen from Eqs. ( H5l) and ( l49l) . Indeed, without the additivity of the 
energy, the entropy cannot be additive as discussed in Sect. |Vl The basic idea is very simple 
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but profound. The entropy being a state function must only depend on state variables. By 
definition, all state variables of a body must be solely determined by what happens within the 
body. If there are other bodies, their influence on the body must not destroy the additivity 
of all extensive state variables. Let us consider the isolated system consisting of the system 
S and the medium E it is in contact with, which is considered in Sect. El in particular, 
consider the energies in Eq. ( H5|) . We will assume no relative motion between E and the 
medium S to make the discussion simple. The energies Eq, E{t) and E{t) then represent 
the energies of the isolated system, the system and the medium, respectively. However, the 
presence of the interaction energy EQ^^\t) in this equation destroys the required property 
that the right side of Eq. f HS]) is the sum of internal energies of the system and the medium; 
the interaction energy is a property of both bodies. Accordingly, the entropy So{Eq) of Sq 
cannot depend only on the energies of the two bodies separately. The difference represented 
by SQ^^\t) in Eq. (!50|) depends on both bodies. Thus, the entropy of the isolated system 
cannot be a sum of entropies, each representing the entropy of one of the bodies alone. 
In other words, generalizing this result to a collection of bodies bj forming a body B, we 
conclude that the entropy ^(B) of B cannot be expressed as a sum of entropies, each term 
S{hj) representing the entropy of the body hj 

5(B) ^^S(b,). 

j 

However, the additivity requires not only that the sum of the energies of various parts b^ of 
a body must yield the energy of the body B itself 

E{B) = Y,E{h,), 
j 

which is satisfied if the linear sizes of the bodies are large compared with the range of 
interaction, but also requires that various bodies are quasi- independent, which is satisfied 
if the linear sizes of the bodies are large compared with thecorrelation length. Unless thess 
properties hold, the entropy cannot be expressed as a sum of the entropies of its parts, with 
each entropy depending only on the properties of the part alone. It is the latter property 
that makes entropy a state function. 
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B. Concept of Internal Equilibrium 



For a body (or its parts, such as subsystems) out of equilibrium, the entropy usually has 
an explicit dependence on time in addition to the implicit time-dependence 

5(B) = 5(X,I,t); 

the latter arises from the dependence of the entropy on state variables X(t) and l{t) that 
explicitly depend on time; we will suppress this time-dependence in the following for nota- 



tional simplicity. (We will assume, as discussed in 33|], that at least one observable, which 
we take to be the number of particles N is held fixed and is not contained in X.) Thus, such 
an entropy will continue to change (increase if the body is isolated) for fixed state variables. 
For a homogeneous body, this will happen if the state variables in X and I do not uniquely 
specify its state. In other words, there may be other internal variables not contained in I. 
This is a connsequence of a simple generalization of Theorem SI The other possibility is 
that the body is not homogeneous. If the isolated body is not in equilibrium, its entropy 
continues to increase according to the second law in Eq. ([2]). However, if it happens that 
the entropy of the body does not exphcitly depend on time 

5ie(B) = ^(X,I), 

then its entropy must be at its maximum for fixed X and I. This is how we had introduced 
the concept of internal equilibrium in I: in this state, the entropy of a system, which is not in 
equilibrium with the medium E(Yo,0), has no explicit time-dependence. It can be brought 
in contact with another medium E', representing the medium S(Yie, Aie), so that it remains 
in equilibrium with the new medium. In other words, there is no difference between a body 
in internal equilibrium with the medium S(Yo,0), and the body in equilibrium with the 
medium E(Yie, Aie); in the latter case, the medium ensures to keep the averages Xie and 
IiE fixed. In both cases, the entropy is maximum for the fixed values of their state variables 
Xie and Iie- Accordingly, all properties of a body in equilibrium also hold for a body in 
internal equilibrium at each instant. For example. Theorem [T] established in Sect. Ill Dl shows 
that there cannot be any relative motion within such a body. The only possible motion is a 
uniform translation and a rigid-body rotation. As a consequence, as discussed in Sect. IIIDl 
there is no viscous dissipation within the body in internal equilibrium. Thus, while there 
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FIG. 5: Schematic form of the entropy S{t) (solid curve) and Seqit) (dotted upper curve) in a 
cooHng experiment. Both curves start from S' at t = and end at S''^ as t ^ oo. The entropy 
difference AS{t) is shown by the leftmost downward pointing arrow and AeS'(t) by the rightmost 
pointing arrow. The irreversible entropy generation AiS{t) > is shown by the small upward 
pointing arrow. 

is no relative motion between the system and S(Yie, Aie), there will in general be relative 
motions between the system and the medium S(Yo,0). This result can be generalized to 
the following statement: The only source of viscous dissipation in a collection of bodies is 
then due to relative motions between different bodies, when each of which is postulated to 
be in internal equilibrium. 

It was assumed in I that the system S at each instant remains in internal equilibrium as it 
goes from one equilibrium macrostate Ai' in contact with a medium S(Yq, 0) at time t = 
to another equilibrium macrostate Ai^ in contact with a medium S(Yo,0) as t — ?■ oo. Let 
us try to understand the behavior of its entropy S{t) during this process, which can then be 
extended to a subsystem in internal equilibrium. The field variables for the two macrostates 
are W = (Y' = Yo,A' = 0) and W" = (Y = Yo,A = 0), respectively. As discussed 
several times, the intermediate macrostate A/'(t) at some instant t can be thought of as an 
equilibrium state of the system E in contact with a medium E(Yie = Y(t), Aje = A(t)) 
after being disconnected from the medium S(Yo,0). In this equilibrium state, the system 
has a well-defined field W(t) = (Y(t), A(t)) and the entropy S{t), which is shown by the 



78 



solid curve in Fig. [51 However, this equilibrium state of the system does not represent the 
equilibrium state of the system in contact with either the medium S(Yo, 0) or the medium 
S(YiE,0). If we bring the system in contact with the medium E(Yie,0), its entropy will 
increase as the internal variables relax towards the new equilibrium, to be denoted by Ai {t) 
with fixed Xie- The equilibrium entropy of the system in Ai(t) will be denoted by 5'eq(t), 
which is shown by the dotted curve in Fig. [5l 

Because of our assumption of internal equilibrium, S{t) represents the maximum possible 
entropy for the macrostate Af{t) at t. The change in the entropy AS{t) is given by 

AS{t) = S{t) - S', 

where S' is the equilibrium entropy in the macrostate Ai' at t = 0, and is shown by the 
left-most arrow in Fig. O The total entropy change in the entire process is 

A5oo = S'- S', 

where as the equilibrium entropy in the macrostate A^" as t — t- oo. 

The entropy difference AS{t) can always be broken into two parts in analogy with the 
two terms in eq. fll26p : 

AS{t) = AeS{t) + AiSit), 

with AcS{t), see the right-most arrow in Fig. [5l representing the entropy exchange with the 
medium S(Yo,0) and A^S^ > the irreversible entropy production within the system, see 
the middle small arrow in Fig. [5l It does not matter whether the process is reversible or 
irreversible, the entropy change is always AS(t). In the former process, AcS(t) = AS{t), 
AiS{t) = and S{t) = Seq{t). In the latter process, 

A^S{t) < AS{t) 

and there is irreversible entropy production during the passage from A4' to Af{t), 

\S{t) > 

in accordance with Eq. (11281) . We thus see that irreversible generation of entropy has 
been accounted for in our approach; see also I and Eq. (11281) . We can easily extend the 
above discussion to the subsystems, where we now also have the possibility of irreversible 



79 



entropy generation through relative motions between them, which is also accounted for in 
our approach. 

We have found (Theorem Hj) that while a body is in internal equilibrium in the presence of 
internal variables, it is not in internal equilibrium when internal variables are not considered 
in its description. In other words, the entropy of a body, see the solid curve for S{t) in Fig. 

El 

5(B) = 5ie(B) = 5(X,I) (136) 

has no explicit time-dependence; however, it becomes an explicit function of time if expressed 
only in terms of observables: 

5(B) = 5(X,t). (137) 

This theorem shows how important the internal variables are in describing glasses. Because 
of the presence of t, this entropy will continue to increase and approach Seq{t), shown in 
Fig. m if we keep X fixed. 

As experimentalists can only control the observables in X, it is the expression of the 
entropy in Eq. f ll37p that is experimentally relevant. Because of the explicit t-dependence, 
the body will continue to relax even if the observables are held fixed by isolating the body 
from its surroundings and bringing it in contact with the medium E(Yie, 0). This relaxation 
occurs because the internal variable I is no longer held fixed. It will continue to change with 
time until finally the entropy reaches its maximum possible value for fixed X. We can apply 
Theorem [3] for an isolated system to our body. The conclusion is that the maximum of the 
entropy occurs when A for the internal variable I is identically zero, which is why we had 
selected this particular medium. The final value of I is some value leq(X). The body in this 
equilibrium state is the fully relaxed body with its entropy given by 

'S'eq(B) = 5'cq(X); 

this entropy is shown as S'eq(t) in Fig. [51 The affinity A of the internal variable is identically 
zero, and the fields in this state are given by Y = Yie- This equilibrium state of the isolated 
body should not be confused with the state of the body in internal equilibrium in the presence 
of the internal variable I, as shown in Eq.( ]136p . This body has I different from Icq(6X), and 
the corresponding affinity A different from zero. 
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C. Thermodynamic Potentials 

The additivity of the entropy and energy allows us to treat our system as a collection 
of quasi-independent subsystems so that we can develop the thermodynamics of an inho- 
mogeneous system. Under the mildest possible assumption that the medium is in internal 
equilibrium so that its field variables Wk are defined via Eq. (!98l) . but the system itself is 
(or the subsystems are) not necessarily in internal equilibrium, we are able to identify the 
thermodynamic potential for the system. The field variables of the system in internal equi- 
librium vary with time and, therefore, they differ from the medium's field variables. They 
become identical only when the system is in equilibrium with the medium. It is obvious 
from Eq. (1^ that the product 

WkZk, no summation implied 

is dimensionless because the entropy is defined as dimensionless in this work. Accordingly, 
the (dimensionless) entropy S{Z) can be expressed in terms of these dimensionless products. 
It follows then that resulting thermodynamic potentials will contain these products, as 
shown in Sect. IX B[ The actual form of the thermodynamic potential function depends 
on the choice of the set C, and requires field variables to maintain the dimensions of each 
term appearing in it. The fields variables in the thermodynamic potentials are the field 
variables of the medium, and not of the system as the latter are not even defined when 
the system in not in internal equilibrium; see for example, Eq. fj56|) or fl63l) . This is true 
whether we consider the system or any of its many subsystems; the latter are considered 
in Sect. IX B[ This is surprising since one knows that, in equilibrium thermodynamics, the 
thermodynamic potentials are state quantities. Therefore, they must contain all quantities 
related to the system. It just happens that in equilibrium the field variables are identical 
to the medium's field variables. Therefore, our thermodynamic potentials reduce to the 
standard thermodynamic potentials in equilibrium. In this case, they become state functions. 
But this is not true when the system is not in equilibrium. What is surprising is that 
this result remains true even when the system or subsystem is in internal variable. The 
thermodynamic potentials always contain fields of the medium and not of the system or 
subsystem. These functions have the required thermodynamic property that they can never 
increase in any spontaneous process, as seen from Eq. fllOOp or (11021) . One can easily see that 
this property is a consequence of the convexity property of the thermodynamic potentials. In 
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contrast, a function obtained from thermodynamic potentials by replacing medium's fields 
by the body's field, when the latter is in internal equilibrium, does not have this required 
thermodynamic property. An example of this is the function ^{t) or G{t) from g(t) or G{t) 
given in Eq. or flS^ . 

While each subsystem is in internal equilibrium, so that there can be no irreversible 
processes inside them, the irreversibility occurs due to relative motions going on among 
them. Because of the internal equilibrium, each subsystem has its own internal temperature 
Tk(t), pressure Pkit) or other conjugate field variables. How are these temperatures and 
pressures related to the internal temperature T{t) and pressure P{t) introduced in I? To 
answer these questions, we consider their definitions in the C frame: 

1 _ dS_ P{f)_ _ dS_ 
1 dsk Pk{t) dsk 



see Eq. (\122\\ for the definition of Tk{t) and Pk{t)- Let us first consider the temperature. 
Introducing 



we have 



Thus, we have 



Introducing 



we find that 



dE 

rk{t) = ^,T.rk{t) = 1, 



dE ^^dejdck ^"mty 
Tit) ~ ^kTuity 



P{t) ^ rl{t)P,{t) 



T(t) Tkit) 
The same exercise can be carried out for other state variables. Introducing 

for the l-th state variable Zi, we find that 

' - T(t) - ^k nit) ■ 
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D. Contrast with Local Non-equilibrium Thermodynamics 



Our approach differs from the traditional local non-equilibrium thermodynamics due to 
de Groot |9l-ll2|. briefly discussed in Sect. IIIH in three important ways. 



1. The first one relates to the principle of additivity of energy. The energy in the local 
non-equilibrium thermodynamics is not the sum of the energies of its various part due 
to the presence of the mutual interaction energies expressed by ilj{r)dV as seen from 
Eq. (jlj). Indeed, usually one applies local equilibrium to a volume element dV, what 
is traditionally called a physically infinitesimal volume in that, while it contains a 
large number of atoms, the corresponding volume is still infinitesimally small. Such a 
volume is conventionally called a "particle" (not to be confused with our usage which 
refers to an atom or molecule). For all practical purposes, it is indeed considered 
as a limit dV 0. It is evident that one must then consider the interaction energy 
■ilj{r)dV to account for the interaction of this volume with the rest of the system. Thus, 
the local non-equilibrium thermodynamics takes the additivity of entropy is taken as 
postulate even when the energy is not additive. 

2. Because the volume element dV is treated infinitesimal in the local theory, all ther- 
modynamic quantities are also treated as continuous in space, while this continuity is 
not a prerequisite in our approach. 

3. Our approach also differs from the local non-equilibrium theory in the form of the 
thermodynamic potentials. The volume element is considered to be in internal equi- 
librium from the start so that it has a well-defined temperature, pressure, etc. These 
fields are used in identifying the thermodynamic potentials. For example, the Gibbs 
free energy of the "particle" is taken to be ^{t)dVeven if the system consisting of such 
"particles" is not in equilibrium with the medium. 

4. Another imporatnt difference between the two approaches is that the reversible entropy 
change and the irreversible entropy generation in each subsystem also depends on the 
equilibrium state of the system. The irreversible entropy generation in the local theory 
depends only on the current local properties. 
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These differences make our approach very different from the local non-equilibrium ther- 
modynamics due to de Groot j9-12|. 
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Appendix A: Relation Between Lab and Body Frames 

Let us consider a particle of mass m in the lab frame C, where it has a velocity V£ and 
the potential energy U. The Lagrangian Lc in this frame is given by 

Lr = -mv'i — U. 
^ 2 ^ 

Let us transform to a rotating frame C which is moving with a velocity V and rotating with 
an angular velocity fl. The velocity of the particle in C is given by and is related to 
according to 

V£ = vc+V + nxrc, (Al) 

where rc is the coordinate of the particle in the C frame, and is related to the coordinate rc 
of the particle in the lab frame C by 

rc = R + rc. 

In the following, it would also be convenient to consider a nonrotating frame X, which is 
only moving with the velocity V with respect to the lab frame, but whose origin coincides 
with that of C. The Lagrangian in the frame C is given by 

1 9 

Lc = -m(vc + V + X r^Y - U 

= ^m\l + ^m{n X r^)^ + mv^ • V+mvc ■ fl x rc+mV ■ fl x rc—U, (A2) 

in which we have omitted |mV^, which is a total time derivative. 
The canonical momentum is obtained as 

dLc 

Pc = — = m(vc+V + ri X rc), 

OVc 
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so that the energy of the particle becomes 

Ec = pc-vc- Lc = ^mvc - x r^)^ - mV ■ Q x rc+U. (A3) 

Expressing in terms of V£ from Eq. flAl|) . we find that 

Ec = Ec + -mV^ - P£ • V - mrc x V£ ■ =Ec + -mV^ - p^-V - m ■ f2, (A4) 
2 2 



where 



-E'z: = ^"^v^ + f/, p£ = mv£, m =mrc x V£. (A5) 



Incidentally, we also note that 

Pc= Pc- 

Thus, 

m = rc X p£ = rc X pc, 
which explains why there is no need to use £ or C as a subscript in m. 

Appendix B: A Rotating and Translating System 

Let us now extend the previous calculation for a single particle to a system of particles 
of total mass M at a given instant t. The system is also characterized by the number of 
particles and its volume. For specificity, we use N and V{t) to denote these quantities. 
The notation should not mean that we are only considering the system S here. The system 
we have in mind is any generic system. We assume that this system is translating with 
a velocity V(t) and rotating with an angular velocity Q,{t) as a whole, both of which can 
change in time. For notational simplicity, we will suppress the explicit t-dependence of 
various quantities here. We focus on one particular instant t. We take the center of mass 
of this system to coincide with the origin of C, so that C is fixed to the body and rotating 
with it. In this case, C represents the center of mass frame. We now sum Eq. f lA4|) over all 
particles, with the result that it is replaced by 

Ec = Ec + ^-Pc-y-M-n, (Bl) 

where we have introduced 

P=MV, Pc = Y,^^c, M = J]rc X p^ = J]rc X PC. (B2) 
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Here, Mc is the intrinsic angular momentum of the system of particles in the C frame. The 
summation in the above formulas, which are also applicable to non-uniform rotation and 
translation of the center of mass frame C, is a sum over all the particles in the system. The 
value of Eq in the C frame represents the internal ener gy of the system. The equation fIBip 



generalizes the result given by Landau and Lifshitz [26|, see their Eq. (39.13)] to the case 
when V 7^ 0. The present generalization is not limited to constant rotation and translation. 
We can express M as follows: 

M =^mrc X (vc+V + f2 x rc) = ^mrc x vc + (y^mrg) x V+^mrc x (ft x vq), 

in which the second sum on the right vanishes for C, the center of mass frame. For the same 
reason, the third term in the second equation in Eq. (1A3P does not contribute to the energy 
of the system. Thus, we find 

M=y~^mrc X Vc + ^mrc x (ft x rc) (B3a) 

Ec = IY^^^^c - l^nifl X r,f+U. (B3b) 

We see that 

M ■ n = J^mrc ■ (vc x fl) + ^m(ri x rc f. (B4) 

The first term vanishes when the motion is radial. Thus, it is the contribution to the energy 
from the relative transverse motion in the C frame and will vanish if the system is stationary 
in the this frame. The latter condition is met when the system is in internal equilibrium in 
accordance with Theorem [H 

For the center of mass frame C, P = P£, so that the energy of the system in the frame C 
is given by 

p2 

Ec = Ec M n, (B5) 

2M ^ ^ 



which can be rewritten as 



Ec = Ec + ^ + M-n (B6) 



Using Eq. flB3bl) . we find that 

>2 

2M ' 2 



Ec = Ec + ^ + ^5^m(ll X rc)' + J^mrc ■ (vc x n) , (B7) 
in which the last term is the contribution of the transverse motion. 
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The energy of the system in the X frame, in which the system has only rotation, is given 



by 



Et = Ef 



2M 



(B8) 



It is obvious that Ex does not depend on the velocity V. Thus, 

Ec = Ex -M-n 

does not depend explicitly on V. We can take Ec as a function of Ex and fl, from which it 
follows that 

'^''') =0, (B9) 

/ Ex,V,N,n 

= -M. (BIO) 



dV J 

dEc 



on 



Exy,N 

All the above results are for a particular microstate of the system undergoing a Hamiltonian 
dynamics. To obtain thermodynamics of the system, we need to average the above equations 
over all microstates using their probabilities, which will be taken up in Sect. IVIIII The 
averaging takes into account the stochastic nature of the dynamics, which has not been 
considered in both appendices. 
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